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Introduction 


Description  of  Harbor 

Barbers  Point  Harbor,  Hawaii,  is  located  on  the  southwest  coastline  of  the 
island  of  Oahu  (Figure  1).  Figure  2  shows  the  present  harbor  complex  (subse¬ 
quently  referred  to  as  Plan  la),  consisting  of  an  entrance  channel,  deep-draft 
harbor,  barge  basin,  and  a  private  resort  marina  (often  referred  to  as  the  West 
Beach  Marina).  The  entrance  channel  is  140  m  (450  ft)  wide,  945  m  (3,100  ft) 
long,  and  13  m  (42  ft)  deep,  referenced  to  mean  lower  low  water  (mllw).  The 
deep-draft  harbor  basin  is  12  m  (38  ft)  deep,  700  m  (2,300  ft)  wide,  and  640  m 
(2,100  ft)  long,  covering  an  area  of  0.37  sq  km  (92  acres).  Rubble-mound  wave 
absorbers  line  approximately  1,128  m  (3,700  ft)  of  the  inner  shoreline  of  the 
harbor  basin.  The  barge  basin,  located  just  seaward  of  the  harbor  on  the  south 
side  of  the  entrance  channel,  has  little  shelter  from  incident  wave  energy.  It  is 
67  m  (220  ft)  wide,  396  m  (1,300  ft)  long,  and  7  m  (23  ft)  deep.  West  Beach 
Marina  was  built  on  the  west  side  of  the  deep-draft  harbor.  It  shares  the  same 
entrance  channel  with  the  harbor,  is  4.6  m  (15  ft)  deep,  and  covers  approximately 
0.08  sq  km  (20  acres).  The  marina  was  designed  to  accommodate  350  to 
500  pleasure  boats. 


Figure  1 .  Project  location 
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Figure  2.  Aerial  photograph  of  Barbers  Point  Harbor,  Oahu,  Hawaii 


Background 

Construction  of  Barbers  Point  Harbor  was  initiated  in  March  1982  and 
completed  in  December  1982  at  a  total  cost  of  $47.2  million.  Six  years  later,  on 
behalf  of  the  State  of  Hawaii,  the  U.S.  Army  Engineer  District,  Honolulu, 
authorized  an  assessment  of  the  engineering  feasibility,  costs,  and  physical 
requirements  for  potential  future  modifications  of  Barbers  Point  Harbor  to 
service  vessels  larger  than  the  design  vessel  for  which  the  harbor  was  originally 
constructed.  Barbers  Point  Harbor  was  originally  designed  to  accommodate  a 
Matson  Navigation  Line  (MNL)  Enterprise-class  general-cargo  vessel  with  the 
following  dimensions: 

Length  =  219  m  (720  ft) 

Beam  =  29  m  (95  ft) 

Loaded  draft  =  10.4  m  (34  ft) 

General-cargo  vessels  dominated  shipping  in  Hawaii  in  the  1960s;  however, 
offloading  of  such  vessels  was  very  inefficient.  Competition  caused  companies 
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to  move  towards  larger  vessels  that  can  be  more  quickly  offloaded,  such  as 
containerships  and  roll-on/roll-off  vessels  (e.g.,  car-carriers).1 

Current  and  potential  future  users  of  the  harbor  would  like  the  capability  for 
fully  loaded  vessels  to  enter  and  leave  the  harbor,  but  safe  vessel  draft  is 
restricted  by  existing  harbor  and  entrance-channel  depths.  Numerical 
simulations  of  Barbers  Point  Harbor  were  conducted  in  1988  using  an  American 
President  Line  (APL)  C8-class  containership  with  the  following  dimensions: 

Length  =  240  m  (787  ft) 

Beam  =  30.5  m  (100  ft) 

Loaded  draft  =  10.7  m  (35  ft) 

That  study  indicated  a  C8-class  containership  could  safely  enter  the  harbor  under 
ideal  conditions  (daytime,  no  waves,  no  current)  with  tug  assistance. 

The  U.S.  Army  Engineer  Division,  Pacific  Ocean  (POD),  and  the  Hawaii 
Department  of  Transportation  (HDOT),  Harbors  Division,  sponsored  a 
navigation  study  of  Barbers  Point  Harbor  at  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  between  1990  and  1992.2  Numerical  and  physical 
model  studies  of  harbor  response  to  modifications  to  the  entrance  channel  and 
harbor  geometry  were  conducted.  Physical  model  navigation  experiments  with 
an  APL  C9-class  containership  were  also  conducted.  The  study  recommenda¬ 
tions  included  deepening  the  harbor  from  11.6  m  (38  ft)  to  13.7  m  (45  ft), 
deepening  the  entrance  channel  from  12.8  m  (42  ft)  to  14.9  m  (49  ft),  and 
widening  the  outer  305  m  (1,000  ft)  of  the  entrance  channel  from  137  m  (450  ft) 
to  228  m  (750  ft). 

Table  1  shows  the  increasing  size  of  the  design  vessels  used  in  previous 
studies  and  the  size  of  the  design  vessels  used  in  the  present  model  study.  The 
length  of  the  design  vessel  has  increased  approximately  20  percent  while  the 
draft  has  increased  over  30  percent.  It  is  more  economical  to  the  shipping 
companies  if  fully  loaded  vessels  enter  and  leave  the  harbor.  However,  for  a 
vessel  to  enter  and  leave  the  harbor,  there  must  be  a  safe  distance  between  the 
bottom  of  the  vessel  keel  and  the  ocean  bottom  (underkeel  clearance).  The  fully 
loaded  design  vessels  used  in  this  present  study  could  not  enter  and  leave  the 
now-existing  harbor  configuration  safely.  The  components  of  underkeel 
clearance  are  shown  in  Figure  3.  Although  all  components  of  underkeel 
clearance  are  shown,  only  wave  response  and  squat  allowance  were  addressed  in 
this  study.  Changes  in  water  density  would  not  be  applicable  to  this  project. 

Differences  exist  between  the  State  HDOT  and  the  Corps  POD  underkeel 
clearance  criteria.  Presently,  the  State  HDOT  is  using  a  0.6-m  (2-ft)  underkeel 


1  Permanent  International  Association  of  Navigation  Congresses  (PIANC).  (1987). 
“Development  of  Modem  Marine  Terminal,”  Supplement  to  Bulletin  No.  56,  Brussels,  Belgium, 

2  Briggs,  M.,  Lillycrop,  L.,  Harkins,  G.,  Thompson,  E.,  and  Green,  D.  (1994).  “Physical  and 
numerical  model  studies  of  Barbers  Point  Harbor,”  Technical  Report  CERC-94-14,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 
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Table  1 

Summary  of  Design  Vessels  Used  in  Previous  and  Present  Studies 

Description 

Length,  m  (ft) 

Beam,  m  (ft) 

Fully  Loaded 
Draft,  m  (ft) 

Date  of  Study 

MNL  Enterprise-class 

219 

29 

10.4 

1963 

vessel 

(720) 

(95) 

(34) 

APL  C8-class 

240 

30 

10.7 

1988 

containership 

(787) 

(100) 

(35) 

APL  C9-class 

262 

32 

11.9 

1990 

containership 

(860) 

(106) 

(39) 

Present  study: 

259 

32 

13.7 

Present  study 

modified  Bunga  Saga 

(850) 

(106) 

(45) 

1997 

Empat  bulk-cargo 
carrier  and  APLC9- 
class  containership 

Figure  3.  Components  of  underkeel  clearance 


clearance  criterion  in  the  harbor  that  translates  to  a  1.8-m  (6-ft)  underkeel 
clearance  in  the  channel.  The  Corps  POD  criterion  is  more  conservative,  calling 
for  a  1.2-m  (4-ft)  and  2.4-m  (8-ft)  underkeel  clearance  in  the  harbor  and  entrance 
channel,  respectively.  The  underkeel  clearance  is  shown  graphically  in  Figure  4 
for  both  the  State  HDOT  and  the  Corps  POD  criteria. 
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Study  Objectives 

Objectives  of  this  study  include,  for  specific  design  vessels,  a  determination 
of  the  optimum  vessel  draft/entrance-channel  depth  combination  that  can  safely 
transit  the  entrance  channel  and  harbor.  To  determine  the  optimum  vessel 
draft/entrance-channel  depth  combination,  two  design  vessels  were  chosen:  (a)  a 
bulk-cargo  vessel  with  a  “box”  shape  (modified  Bunga  Saga  Empat)  and 
(b)  a  retrofitted  sleeker  vessel  that  had  been  used  in  the  1990  study  (APL  09- 
class  containership)  (Table  1).  Simulation  of  these  vessels’  transit  through  the 
entrance  channel  was  made  with  selected  wave  conditions,  and  the  underkeel 
clearance  was  measured.  Results  of  the  underkeel  clearance  criteria  data  for 
determining  the  optimum  vessel  draft/entrance-channel  depth  combinations  are 
presented  in  Appendixes  A  through  F. 

Subsequent  to  the  completion  of  the  physical  model  navigation  study 
described  above  for  determining  the  optimum  vessel  draft/entrance-channel 
depth  combination,  POD  requested  WES  to  conduct  numerical  investigations 
that  would  (a)  analyze  impacts  of  the  north  jetty  on  wave  climate  in  the  barge 
basin  and  (b)  document  harbor-oscillation  characteristics  for  entrance  channel 
and  harbor  depths  of  13.4  m  (44  ft)  and  12.8  m  (42  ft),  respectively,  and  14.3  m 
(47  ft)  and  13.7  m  (45  ft),  respectively,  for  a  planned  harbor  expansion  to  183  m 
(600  ft)  wide  and  335  m  (1,100  ft)  long.  Results  of  these  additional 
supplemental  numerical  investigations  to  (a)  analyze  impacts  of  the  north  jetty 
on  wave  climate  in  the  barge  basin  and  (b)  document  harbor-oscillation 
characteristics  are  presented  in  Appendixes  G  and  H,  respectively. 


Figure  4.  State  HDOT  and  Corps  POD  underkeel  clearance  criteria 
shown  graphically 
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2  Physical  Model 


This  chapter  describes  equipment  and  techniques  employed  in  the  physical 
model  to  reproduce  prototype  conditions  at  Barbers  Point  Harbor.  Model  scaling 
laws  were  used  to  ensure  similarity  between  the  prototype  and  model  for  both 
the  harbor  and  required  vessels.  The  following  equipment  was  utilized  to 
reproduce  and  verify  that  the  proper  incident  wave  conditions  were  simulated  in 
the  model:  a  wavemaker  to  generate  selected  incident  waves  and  wave  gauges  to 
ensure  that  the  proper  wave  heights  were  being  generated.  A  high-resolution 
video-tracking  system  was  used  to  record  vessel  motions  during  transits  both  into 
and  out  of  the  harbor. 


Model  Design 

An  undistorted,  three-dimensional  model  of  Barbers  Point  Harbor  (Figure  5) 
was  constructed  in  1991  at  a  model  to  prototype  scale  Lr  =  1:75,  in  accordance 
with  well-known  Froude  scaling  laws.  The  nearshore  area  extends  to  the  30.5-m 
(100-ft)  mllw  contour  and  includes  approximately  1  km  (3,500  ft)  on  either  side 
of  the  entrance  channel.  Total  area  of  the  model  is  over  1,000  sq  m  (11,000  ft2). 
This  model  scale  was  selected  to  allow  proper  reproduction  of  significant  harbor 
features.  Model  and  prototype  lengths  scale  as  L„  areas  scale  as  L2  and  time 
and  velocity  scale  as  \/Lr  Details  of  the  model  design  can  be  found  in  Briggs 
et  al.  (1994).1 


Model  Appurtenances 

Wavemaker 

At  slow  vessel  speeds,  the  principal  component  of  underkeel  clearance  is  the 
wave-induced  vessel  response.  To  simulate  realistic  wave  conditions  found  at 
Barbers  Point,  waves  were  generated  with  a  plunger  wavemaker  that  is  capable 
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2  Briggs,  M.,  Lillycrop,  L.,  Harkins,  G.,  Thompson,  E.,  and  Green,  D.  (1994).  “Physical  and 
numerical  model  studies  of  Barbers  Point  Harbor,”  Technical  Report  CERC-94-14,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 
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Figure  5.  Barbers  Point  Harbor  physical  model 


of  generating  irregular,  unidirectional  wave  conditions.  It  is  called  a  “plunger” 
wave  machine  because  the  wave  board  is  constructed  in  the  shape  of  a  trapezoid 
and  is  plunged  into  the  water  column.  The  shape  of  the  wave  paddle  is  shown  in 
Figure  6. 

Simulating  the  appropriate  wave  conditions  is  a  process  consisting  of  the 
following  steps.  The  proper  wave  parameters  for  the  model  wave  climate  are 
calculated,  and  a  control  signal  representing  the  wave  paddle’s  displacement 
from  its  mean  position  as  a  function  of  time  is  written  digitally  to  tape.  For  the 
plunger  wave  machine,  this  consists  of  a  single  time  series  since  the  entire  wave 
paddle  moves  in  unison.  The  waves  are  generated  and  analyzed,  and  the  control 
signal  time  series  is  then  modified  to  increase  or  decrease  the  wave  height.  No 
correction  is  generally  needed  for  the  wave  period. 


Wave  gauges 

Gauges  were  placed  in  the  model  to  measure  wave  conditions  generated  by 
the  wavemaker  so  that  needed  adjustments  to  the  wave  signal  could  be 
calculated.  Three  parallel-wire  capacitance  gauges  were  used  to  measure  model 
wave  heights  and  were  located  along  the  9-m  (30-ft)  depth  contour  (Figure  7). 
To  ensure  accurate  measurements  of  the  water-surface  elevation,  wave  gauges 
were  calibrated  each  day. 
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Figure  6.  Plunger-type  wave  generator  used  in  physical  model  study 

Computer  support 

Four  computer  systems  were  used  throughout  the  study.  The  tasks  that  these 
computers  accomplished  include  generation  of  the  control  signals  for  simulating 
the  wave  fields,  collecting  wave  data,  tracking  the  design  vessel,  and  analyzing 
the  data. 


Water-level  controller  and  point  gauge 

Water  depth  was  maintained  within  ±0.01  cm  (0.002  ft)  of  the  desired  level 
by  an  automatic  water-level  float  and  an  on/off  solenoid  control  valve.  When  the 
float  fell  below  a  predetermined  position,  water  was  allowed  to  flow  into  the 
model  until  the  appropriate  level  was  again  reached. 


Design  Vessels 

Two  design  vessels,  approximately  the  same  size  but  with  different  hull 
shapes,  were  utilized  during  these  present  navigation  simulation  experiments. 
The  primary  design  vessel  used  during  this  study  had  a  bulk-cargo  carrier  hull 
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Figure  7.  Physical  model  layout,  including  wave  generator  and  gauge  locations 


shape  (modified  Bunga  Saga  Empat).  Barbers  Point  Harbor  was  originally 
designed  to  assist  the  island  sugar  industry  by  exporting  bulk  products  to  the 
U.S.  mainland.  Coal  is  also  shipped  through  Barbers  Point.  Through  economic 
analysis,  POD  determined  that  the  largest  savings  in  shipping  costs  could  be 
realized  by  coal  bulk-cargo  transportation.  The  secondary  design  vessel  used  in 
this  present  study  was  a  more  streamlined  APL  C9-class  containership. 


Chapter  2  Physical  Model 


9 


Modified  Bunga  Saga  Empat  bulk-cargo  carrier 

A  model  of  a  bulk-cargo  carrier  (a  modified  Bunga  Saga  Empat)  was 
constructed  to  use  as  one  of  two  design  vessels  (Figure  8).  The  dimensions  of 
the  actual  prototype  vessel  are  as  follows: 

Length  =  229  m  (750  ft) 

Beam  =  32  m  (106  ft) 

Loaded  draft  =  12.2-13.7  m  (40-45  ft) 


Figure  8.  Bulk-cargo  carrier  Bunga  Saga  Empat  from  which  the  primary  model 
was  developed 


However,  POD  decided  to  extend  the  length  of  the  Bunga  Saga  Empat  by  31  m 
(100  ft)  so  that  this  design  vessel  had  an  overall  prototype  length  of  259  m 
(850  ft). 

The  model  of  the  Bunga  Saga  Empat  was  constructed  at  the  same  1:75  scale 
as  the  Barbers  Point  Harbor  physical  model,  from  a  single  blueprint  and  one 
photograph  of  the  prototype  (Figure  8).  The  dimensions  of  the  design  vessel  are 
given  in  Table  2,  and  a  photograph  of  the  model  vessel  is  given  in  Figure  9. 

A  12-volt  battery  powered  the  variable  speed  motor  and  propeller  assembly. 
Motor  speed,  direction,  and  rudder  control  were  all  remotely  operated  using  a 
custom-built  controller.  Vessel  speed  was  accurately  controlled  using  a  digital 
tachometer.  Prior  to  simulating  wave  conditions,  speed  runs  were  made  in  deep 
water.  The  tachometer  setting  that  reproduced  the  selected  speed  was 
recorded  and  used  on  subsequent  runs  when  waves  were  present.  Before  each 
model  run,  the  motor  speed  was  adjusted  to  the  correct  settings  while  the  model 
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Table  2 

Dimensions  of  Primary  Design  Vessel,  Modified  Bunga  Saga 

Empat 

Prototype  Dimensions,  m  (ft) 

Model  Dimensions,  m  (ft) 

Length 

259  (850) 

3.45  (11.3) 

Beam 

32  (106) 

0.43  (1.4) 

Loaded  draft 

13.7  (45) 

0.18(0.6) 

Figure  9.  Physical  model  of  the  primary  design  vessel,  the  modified  Bunga 
Saga  Empat  bulk-cargo  carrier 

ship  was  held  stationary.  The  vessel  was  then  released  at  the  proper  time. 
Differences  in  vessel  speed  are  attributed  to  wave  conditions  and  proximity  to 
the  bottom. 


APL  C9-class  containership 

A  model  of  an  APL  C9-class  containership  had  been  constructed  for  previous 
physical  model  experiments  conducted  at  WES  in  1990,  and  that  vessel  was 
retrofitted  as  the  second  design  vessel  for  use  in  this  present  study.  Unlike  the 
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modified  Bunga  Saga  Empat  model,  the  APL  C9-class  containership  model  has  a 
much  sleeker  hull  shape  (Figure  10).  For  navigation  experiments,  the  containers 
stacked  on  the  deck  were  replaced  with  equivalent  flush-mounted  weights  and 
reflective  video-tracking  balls. 


Figure  10.  Physical  model  of  the  secondary  design  vessel,  the  APL  C9-class 
containership 


Model  Vessel  Calibration 

Modified  Bunga  Saga  Empat  bulk-cargo  carrier 

The  primary  design  vessel  was  the  modified  Bunga  Saga  Empat.  To  satisfy 
true  Froude  scaling  laws,  a  model  must  be  geometrically  (both  in  shape  and 
size),  kinematically,  and  dynamically  similar  to  the  prototype.  Since  detailed 
specifications  of  the  vessel  could  not  be  obtained  from  Hyundai  Heavy  Industries 
(shipbuilder  of  the  Bunga  Saga  Empat),  a  numerical  model  was  utilized  to  calcu¬ 
late  the  dynamic  and  kinematic  vessel  characteristics. 

The  first  step  in  calibrating  the  vessel  was  to  calculate  the  weights  required 
for  each  draft.  To  calculate  the  weights,  the  shape  of  the  hull  had  to  be  repro¬ 
duced  in  the  numerical  model  to  determine  the  weight  of  the  displaced  water  for 
the  different  draft  conditions.  The  total  weight  needed  to  reach  each  draft  is 
equal  to  the  weight  of  water  displaced.  The  vessel  should  have  neither  trim 
(difference  in  elevation  between  the  front  and  the  back  of  the  vessel)  nor  heel 
(difference  in  elevation  between  the  sides  of  the  vessel). 
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Ships  have  a  natural  period  of  response  to  waves  that  must  be  reproduced  in 
the  physical  model.  An  unrestrained  vessel  is  free  to  move  in  six  degrees  of 
freedom  (defined  in  Figure  11),  including  three  translational  motions  and  three 
rotational  motions.  For  underkeel  clearance  calculations,  vertical  motions  con¬ 
sisting  of  heave,  pitch,  and  roll  are  exceedingly  important.  Translational  heave 
motion  is  a  function  of  the  vessel  shape  and  the  vertical  distribution  of  the 
weights.  Heave  occurs  when  the  vessel  moves  up  and  down  as  one  unit.  Roll 
and  pitch  motions  are  a  function  of  the  horizontal  and  vertical  placement  of  the 
weights.  Roll  occurs  when  the  vessel  rotates  about  a  line  running  through  the 
center  of  mass  of  the  vessel  along  the  long  horizontal  axis  of  the  vessel.  The 
center  of  mass  is  a  point  at  which  the  entire  weight  of  the  vessel  could  be  placed 
and  the  vessel  would  behave  the  same.  Pitch  motion  is  rotation  through  the 
center  of  mass  along  a  line  through  the  short  horizontal  axis  of  the  vessel. 
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Figure  11.  Six  degrees  of  freedom  for  vessel  subject  to  wave  motion 

The  natural  response  of  an  unrestrained  vessel  can  be  displayed  graphically. 

A  numerical  model  generates  the  response  amplitude  operator  (RAO)  curves  for 
the  particular  vessel  draft  and  water  depth  (Figure  12  ).  The  RAO  curves 
indicate  the  motion  (either  translational  or  rotational)  of  a  stationary  vessel  to 
waves  of  different  wave  periods.  (The  direction  of  the  waves  varies  based  upon 
the  motion,  but  the  directions  are  chosen  to  produce  maximum  translation  or 
rotation.)  The  x-axis  indicates  the  prototype  period  of  vessel  response,  while  the 
y-axis  has  dual  units.  The  two  rotational  motions  (roll  and  pitch)  have  units  of 
degrees/feet,  while  the  translational  motion  (heave)  has  units  of  feet/feet.  The 
values  on  the  y-axis  indicate  the  degrees  of  rotation  for  a  unit  wave  height  for  the 
rotational  motion  or  the  distance  of  vertical  excursion  for  the  translational  or 
heave  motion. 

To  simulate  the  prototype  weight  distribution,  2.26  kg  (5  lb)  of  lead  weights 
were  distributed  throughout  the  model  hull  (Figure  13).  First,  the  vessel  draft 
was  checked  by  mounting  “J-shaped”  point  gauges  to  the  model  vessel 
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Bulk  Carrier,  LqA  =  850  ft,  0  =  0° 


Draft  =  42  ft,  Depth  =  45  ft 


Figure  1 2.  Response  amplitude  operators  for  1 2.8-m  (42-ft)  vessel  draft  in 
13.7-m  (45-ft)  water  depth 

(Figure  14).  The  “J-shaped”  point  gauges  were  mounted  to  the  front  and  rear 
along  the  center  line  of  the  vessel  and  to  the  left  and  right  sides  along  the  hull  at 
midship  position  (Figure  15).  The  weight  of  the  point  gauges  was  accounted  for 
while  the  boat  was  being  drafted.  After  the  correct  drafts  were  reached,  the 
dynamic  response  of  the  vessel  was  simulated  by  displacing  the  vessel  in  the 
appropriate  direction.  The  response  period  of  the  vessel  was  checked  by  using 
stop  watches  and  by  counting  the  number  of  cycles. 


APL  C9-class  containership 

The  secondary  design  vessel,  the  APL  C9-class  containership,  had  been  used 
in  earlier  tests  and  was  refurbished  and  instrumented  for  the  present  study. 
Unlike  the  modified  Bunga  Saga  Empat  that  was  numerically  modeled  to 
determine  the  six  degrees  of  freedom,  the  characteristics  of  the  APL  C9-class 
containership  were  provided  by  the  manufacturer  and  were  reproduced  during 
the  1992  Barbers  Point  Harbor  study.  This  vessel  had  previously  been 
appropriately  ballasted. 
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Figure  1 3.  Distribution  of  lead  weights  throughout  the  hull  of  the  model  of  the 
modified  Bunga  Saga  Empat  bulk-cargo  carrier 
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Figure  14.  “J-shaped”  point  gauge 
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Figure  1 5.  Location  of  “J-shaped”  point  gauges  on  the  model  of  modified  Bunga  Saga 
Empat  bulk-cargo  carrier 
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3  Vessel  Motion  Analysis 


A  software  and  hardware  system  developed  by  Motion  Analysis  Corporation 
was  used  to  collect  and  analyze  the  model  ship  motion.  This  is  a  video-based 
system  that  uses  cameras  and  strobes  to  track  reflecting  balls  in  the  field  of  view. 
By  using  six  Styrofoam  reflecting  balls  on  the  vessel  and  four  Styrofoam 
reflecting  balls  located  around  the  channel,  the  distance  as  a  function  of  time 
between  the  vessel  and  the  bottom  could  be  defined. 

The  0.6-m  (2-ft)  prototype  difference  between  the  State  HDOT  and  the  Corps 
POD  harbor  underkeel  clearance  criteria  translates  to  0.8  cm  (0.32  in.)  at  model 
scale.  The  model  arrangement  had  to  measure  at  least  0.3-m  (1-ft)  prototype 
underkeel  clearances,  which  translates  to  0.4  cm  (0.16  in.)  at  model  scale. 


Motion  Analysis  Equipment  Operation 
and  Model  Vessel  Arrangement 

The  motion  analysis  equipment  uses  Charged  Coupled  Device  digital  cameras 
to  record  the  motion  of  reflecting  targets.  Six  Styrofoam  balls  were  wrapped 
with  highly  reflecting  tape  and  mounted  above  the  ship  hull.  Styrofoam  was 
used  to  keep  the  weight  down  so  that  the  reflecting  targets  did  not  modify  the 
vessel  motion.  The  balls  were  mounted  above  the  ship  hull  using  threaded  steel 
rods  anchored  with  wooden  base  plates.  The  vessel  superstructure  was  removed 
during  data-collection  efforts  to  minimize  masking  of  the  targets  from  the 
cameras.  However,  the  weight  of  the  superstructure  was  accounted  for  with  lead 
weights. 

Three  noncolinear  points  in  space  are  needed  to  define  a  geometric  plane.  If 
these  three  points  move  as  a  rigid  body,  then  the  motion  of  the  plane  can  be 
determined.  Three  balls  were  needed  on  the  vessel  to  determine  the  motion  of 
the  vessel.  However,  if  only  three  balls  were  used,  when  one  of  the  balls  left  the 
field  of  view,  motion  of  the  vessel  could  not  be  determined.  By  using  six  balls, 
at  least  three  balls  stayed  in  the  field  of  view  for  a  longer  period  of  time. 

The  balls  were  located  at  two  elevations  above  the  deck  so  that  the  points  are 
not  linear,  which  makes  the  calculation  of  the  roll  component  possible.  The 
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distance  between  the  upper  ball  and  the  lower  ball  was  determined  from  the 
expected  instrumentation  error.  The  horizontal  separation  between  the  balls 
needed  to  be  at  least  15  cm  (6  in.)  to  give  at  most  0.3  m  (1  ft)  error  in  the 
prototype  measurement  for  a  worst-case  situation  in  which  the  error  in  the 
location  of  the  centroid  of  the  balls  was  off  by  0.25  cm  (0.1  in.)  in  opposite 
directions  at  a  particular  instance  in  time.  The  balls  were  placed  23  cm  (9  in.) 
apart  vertically,  thus  providing  a  worst-case  instrumentation  error  of  0.2  m 
(0.7  ft)  prototype  in  the  rotational  motion.  The  ball  mount  had  to  be  rigid  and 
light;  thus,  there  was  a  tradeoff  between  increasing  the  difference  between  the 
elevation  of  the  balls  to  decrease  the  effect  of  instrumentation  error  and  keeping 
the  mounts  as  stiff  and  light  as  possible.  Schematics  of  the  ball  locations  on  the 
modified  Bunga  Saga  Empat  bulk-cargo  carrier  and  the  APL  C9-class 
containership  are  shown  in  Figures  16  and  17,  respectively.  The  balls  were 
labeled  sequentially  starting  with  ball  number  1  at  the  front  of  the  boat. 


Figure  16.  Location  of  reflecting  balls  on  the  modified  Bunga  Saga  Empat 
bulk-cargo  carrier 


Figure  1 7.  Location  of  reflecting  balls  on  the  APL  C9-class  containership 


Motion  Analysis  Equipment  Positioning 

Before  laboratory  navigation  runs  could  commence,  the  accuracy  of  the 
laboratory  arrangement  of  the  motion  analysis  system  had  to  be  defined.  The 
initial  laboratory  arrangement  of  the  motion  analysis  cameras  used  three  cameras 
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located  approximately  20  ft1  from  the  entrance  channel  center  line  (Figure  18). 
The  six-camera  arrangement  covered  the  entire  length  of  the  entrance  channel. 
Although  the  manufacturer  claimed  accuracies  of  1  mm  (0.04  in.)  over  the  entire 
field  of  view,  accuracies  measured  by  calculating  the  distances  between  rigidly 
mounted  balls  over  the  field  of  view  were  roughly  twice  the  manufacturer  claims 
(at  best).  To  obtain  the  level  of  accuracy  desired,  the  initial  laboratory 
arrangement  of  three  cameras  covering  a  20-ft  field  of  view  was  abandoned 
when  the  accuracy  of  the  system  could  not  be  resolved  to  better  than 
approximately  0.25  m  (1.6  ft)  prototype  overall.  To  obtain  the  level  of  accuracy 
required  by  the  study,  six  cameras  were  used  with  a  combined  field  of  view  of 
15  ft  along  the  entrance  channel. 


Figure  1 8.  Initial  arrangement  of  motion  analysis  cameras 


1  To  convert  feet  to  meters,  multiply  by  0.3048. 
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Early  in  the  study,  it  was  observed  that  the  transition  between  the  harbor  and 
entrance-channel  depths  was  an  area  of  possible  groundings.  This  is  the  area 
that  was  covered  in  detail  with  the  15-ft  field  of  view  (Figure  19).  Since  the 
channel  is  quite  long,  a  second  concise  sampling  effort  was  undertaken  in  the 
channel.  The  camera  arrangement  for  this  effort  is  shown  in  Figure  20. 
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Figure  1 9.  Camera  arrangement  for  the  transition  area  data  collection 
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Figure  20.  Camera  arrangement  for  channel  area  data  collection 


Calibration  of  Field  of  View 

The  cameras  were  placed  in  approximately  the  shape  of  a  circle  positioned  at 
a  6-m  (20-ft)  radius  around  the  center  of  the  field  of  view.  This  radial  distance 
from  the  center  of  the  field  of  view  was  chosen  so  that  each  camera  could  see  its 
entire  portion  of  the  channel.  The  locations  around  the  perimeter  of  the  circle 
were  chosen  so  that  each  camera  could  see  all  six  reflecting  balls  on  the  vessel 
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and  were  not  in  line  with  the  direction  of  vessel  travel  (i.e.,  no  cameras  were 
mounted  at  the  ends  of  the  channel).  After  the  camera  positions  were 
established,  a  calibration  procedure  was  performed  to  obtain  a  three-dimensional 
(3-D)  mapping  of  the  field  of  view.  A  calibration  cube  was  first  placed  in  the 
field  of  view  and  oriented  so  that  all  six  cameras  had  a  clear  view  of  each  of  the 
reflecting  sources  located  at  the  eight  comers  of  the  cube.  Data  were  collected 
and  digitized  for  2  sec  at  60  Hz  for  each  of  the  six  cameras.  A  file  was  then 
created  that  took  into  account  the  six  different  camera  lens’  characteristics  and 
their  associated  distortions  and  the  location  of  the  cameras  to  reproduce  an 
accurate  map  of  the  calibration  cube  position. 

The  calibration  cube  was  0.77  m  (2.52  ft)  long,  0.53  m  (1.7  ft)  wide,  and 
0.62  m  (2.03  ft)  tall.  To  increase  the  accuracy  of  the  mapping  of  areas  beyond 
the  calibration  cube,  a  calibration  wand  was  moved  throughout  the  field  of  view. 
The  calibration  wand  was  1.5  m  (5  ft)  long  with  three  reflecting  balls  located 
0.3  m  (1  ft)  and  1.2  m  (4  ft)  apart. 

To  test  the  accuracy  of  the  final  3-D  mapping,  the  modified  Bunga  Saga 
Empat  was  propelled  through  the  field  of  view.  The  distances  between  the 
reflecting  balls  were  then  calculated,  with  the  first  distance  being  subtracted 
from  the  remaining  values.  The  subtraction  of  the  first  value  was  done  to 
remove  the  actual  distance  so  that  values  relative  to  zero  could  be  evaluated.  For 
a  perfect  mapping  of  the  3-D  space  with  no  instrument  error,  the  values  would  all 
be  zero.  The  results  from  the  distance  measurements  between  balls  1  and  2, 
between  balls  2  and  3,  and  between  balls  3  and  4  are  shown  in  Figure  21.  Here 
are  shown  unedited  data  with  no  quality  checks.  Some  outliers  or  points  far 
from  zero  can  be  seen  that  would  be  removed  if  the  underkeel  clearance  values 
were  to  be  calculated.  The  majority  of  the  points  are  within  ±0.1  in.  There  is  a 
slight  trend  to  the  data  that  on  subsequent  analysis  of  the  underkeel  clearance 
results  was  removed  from  the  data  and  will  be  discussed  later. 


Calculation  of  Underkeel  Clearance 

Virtual  markers 

The  reflecting  balls  used  to  calculate  the  motion  of  the  vessel  were  mounted 
above  water,  and  thus  a  method  to  calculate  the  underkeel  clearance  of  the  vessel 
was  needed.  The  Motion  Analysis  Corporation  software  package  provides  a 
means  to  identify  virtual  markers  at  any  position  on  the  vessel.  Assuming  a  rigid 
body  and  given  the  location  of  three  of  the  balls  on  the  ship,  any  point  along  the 
vessel  keel  could  be  identified.  Although  outlining  the  entire  hull  of  the  ship 
with  virtual  markers  was  possible,  the  amount  of  data  then  collected  would  be 
unmanageable  and  exceedingly  intensive  computationally  and  was  unneeded  to 
accomplish  this  study  purpose. 

By  identifying  key  locations  along  the  hull,  the  maximum  vertical  excursion 
of  the  vessel  can  be  ascertained.  For  the  modified  Bunga  Saga  Empat  bulk-cargo 
carrier,  placing  virtual  markers  as  shown  in  Figure  22  provided  the  maximum 
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Figure  22.  Location  of  virtual  markers  on  modified  Bunga  Saga  Empat  bulk- 
cargo  carrier  for  computing  underkeel  clearance 

vertical  excursion  for  either  the  translational  or  rotational  motions.  For  pure 
pitch,  the  two  markers  along  the  center  line  of  the  vessel  at  the  front  and  back 
provided  the  location  of  maximum  motion,  while  the  off-center  points  captured 
pitch/roll  combined  motion. 

The  shape  of  the  APL  C9-class  containership  warranted  different  virtual 
marker  locations.  For  the  sleek-hulled  vessel,  virtual  markers  located  at  four 
positions  were  needed  (Figure  23).  Again,  two  virtual  markers  placed  one  at  the 
front  and  one  at  the  back  of  the  vessel  captured  pure  pitch  motion,  while  the  two 
virtual  markers  located  along  the  sides  of  the  vessel  at  midpoint  defined  pure  roll 
motion. 


Figure  23.  Location  of  virtual  markers  on  APL  C9-class  containership  for  computing  underkeel 
clearance 


Establishing  location  of  hull  above  the  bottom 

The  cube  and  wand  calibration  established  an  arbitrary  coordinate  system,  but 
defined  the  spatial  relationship  of  the  cameras.  However,  a  coordinate  system 
relative  to  the  physical  model  was  needed.  The  system  software  provides  the 
ability  to  perform  a  hierarchal  translation  and  rotation  of  the  coordinate  system. 
Using  the  stationary  markers  located  around  the  channel  (Figures  18, 19,  and 
20),  the  arbitrary  coordinate  system  was  transposed  into  a  coordinate  system 
where  the  x-axis  points  toward  the  channel  sidewalls,  y-axis  lies  along  the  center 
of  the  channel  pointing  seaward,  and  the  z-axis  points  up.  The  origin  lies  at  the 
transition  location  for  the  first  camera  arrangement  (Figure  19)  and  in  the  center 
of  the  capture  view  for  the  second  camera  arrangement  (Figure  20). 

To  establish  the  vessel  position  relative  to  the  harbor  bottom,  the  following 
technique  was  employed.  Before  navigation  runs  were  performed  for  each  water 
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depth/draft  configuration,  an  initial  data  collection  of  the  vessel  facing  into  and 
out  of  the  harbor  was  performed.  The  vessel  was  positioned  so  that  the  third  ball 
was  located  over  the  coordinate  origin,  and  data  were  recorded  for  5  sec.  Using 
the  static  underkeel  clearance  and  the  initial  position  of  the  vessel,  the  underkeel 
clearance  during  experimental  runs  could  be  calculated. 
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4  Selection  of  Experimental 
Conditions 


The  estimated  cost  to  dredge  Barbers  Point  Harbor  and  channel  to  -43  ft  mllw 
and  -47  ft  mllw,  respectively,  is  $30  million  (POD  1996).  Decreasing  the  depth 
in  the  harbor  and  entrance  channel  by  1  ft  would  save  approximately  $6  million. 
Although  the  2-ft  underkeel  clearance  difference  between  the  State  HDOT  and 
the  Corps  POD  criteria  is  small,  it  translates  to  over  $12  million.  To  determine 
the  underkeel  clearance  criteria  necessary  for  safe  transit  through  the  channel 
and  into  the  harbor,  two  design  vessels  were  evaluated  with  12  combinations  of 
wave  period  and  height  for  two  incident  wave  directions.  Data  were  collected  on 
both  the  inbound  and  outbound  travel  directions. 


Wave  Climate 

Wave  height  and  period 

Information  about  the  wave  conditions  found  at  Barbers  Point  Harbor  was 
extracted  from  a  data  collection  effort  begun  in  1986. 1  A  Datawell  waverider 
buoy  located  approximately  1.6  km  (1  mile)  offshore  (21°  20.1  N  latitude,  158° 
9.0  W  longitude)  at  200-m  (650-ft)  depth  was  installed  to  provide  wave  height 
and  period  data  for  the  incident  deepwater  waves.  This  buoy  was  located  in 
water  deep  enough  to  minimize  the  bathymetric  effects  on  the  measured  waves. 
The  waverider  buoy  uses  a  vertically  stabilized  accelerometer  to  sense  the 
vertical  component  of  the  buoy’s  motion.  Vertical  acceleration  and 
displacement  data  are  transmitted  up  to  50  km  (30  miles)  from  shore.  Data  were 
collected  from  July  1986  through  March  1990  and  are  shown  in  Table  3. 

The  median  wave  height  from  the  offshore  waverider  buoy  falls  in  the  0.9-m 
(3-ft)  wave-height  bin.  The  median  wave  period  falls  in  the  7-sec  period  band; 
however,  waves  with  periods  around  13  sec  are  also  very  common. 


1  Lillycrop,  L.,  Briggs,  M.,  Harkins,  G.  Boc,  S.,  and  Okihiro,  M.  (1993).  “Barbers  Point  Harbor, 
Oahu,  Hawaii,  Monitoring  Study,”  Technical  Report  CERC-93-18,  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 
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Waves  chosen  for  the  navigation  runs  should  represent  typical  conditions 
found  at  Barbers  Point.  Also,  wave  heights  should  not  exceed  those  under  which 
a  harbor  pilot  would  be  willing  to  bring  a  vessel  into  harbor.  The  vessel 
response  is  a  function  of  the  joint  wave  height  and  period.  However,  smaller 
wave  heights  will  produce  smaller  motions  of  the  vessel.  POD  desired  to 
determine  the  maximum  significant  wave  height  under  which  the  harbor  pilots 
would  consider  bringing  a  design  vessel  into  the  harbor  and  in  what  water  depth 
the  wave  heights  they  had  determined  as  being  important  should  be  measured. 
However,  a  specific  water  depth  and  wave  height  could  not  be  determined. 

Based  upon  guidance  from  POD,  the  wave  conditions  defined  in  Table  4  were 
reproduced  in  the  physical  model  for  these  navigation  studies.  Selected 
experimental  wave  conditions  were  defined  at  the  9.1-m  (30-ft)  contour.  To 
ensure  wave  conditions  were  correctly  reproduced  in  the  physical  model,  wave 
gauges  were  placed  along  the  9.1-m  (30-ft)  depth  contour  (Figure  7). 


Table  3 

Buoy  Percentage  Distribution  of  Deepwater  Significant  Wave  Height  and  Period 

Significant 

Wave  Height, 

H,( ft) 

Peak  Period,  sec 

Sum  of 
Observations 

% 

22+ 

20 

17 

15 

13 

11 

9 

7 

5 

10+ 

1 

2 

3 

0.1 

9 

1 

4 

1 

3 

9 

0.2 

8 

5 

6 

4 

6 

2 

2 

2 

27 

0.7 

7 

4 

3 

20 

8 

9 

6 

9 

1 

60 

1.6 

6 

4 

8 

33 

30 

29 

16 

13 

133 

3.6 

5 

1 

22 

79 

92 

49 

35 

35 

6 

319 

8.6 

4 

11 

62 

173 

237 

176 

145 

208 

37 

1,049 

28.2 

3 

10 

47 

182 

384 

307 

306 

527 

87 

1,850 

49.8 

2 

4 

16 

65 

55 

44 

60 

18 

262 

7.1 

1  1 

4 

4 

0.1 

Sum  of 
Observations 

35 

153 

510 

826 

628 

558 

857 

149 

3,716 

% 

1.0 

4.1 

13.7 

22.2 

16.9 

15.0 

23.1 

4.0 

100.0 

The  waves  were  calibrated  at  the  three  gauge  locations  by  first  generating  and 
storing  the  parameters  listed  in  Table  4.  The  file  was  then  used  to  generate 
waves  at  the  same  time  the  wave-gauge  data  were  being  collected.  The  data 
were  analyzed,  and  corrections  were  made  to  the  control  signal  until  the  target 
wave  height  was  found.  Table  5  shows  the  results  for  the  control  signals. 
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Table  4 

Wave  Conditions  Reproduced  in  the  Physical  Model 

Prototype  Wave 

Height,  m  (ft) 

Wave  Condition 
Number 

Prototype  Wave 
Period,  sec 

Gamma  (width  of  the  spectra 
in  the  frequency  domain) 

0.9  (3.0) 

1 

6.0 

3.3 

2 

10.0 

5.0 

3 

14.0 

7.0 

4 

18.0 

10.0 

1 .5  (5.0) 

5 

6.0 

3.3 

6 

10.0 

5.0 

7 

14.0 

7.0 

8 

18.0 

10.0 

2.1  (7.0) 

9 

6.0 

3.3 

10 

10.0 

5.0 

11 

14.0 

7.0 

12 

18.0 

10.0 

Table  5 

Wave  Period  and  Height  from  Control  Signal  Calibration 

Prototype 

Period,  sec 

Prototype  Target 
Wave  Height,  ft 

Measured  Wave  Height,  ft 

Percent  Error,  % 

6 

3.0 

3.02 

0.75 

10 

3.0 

2.99 

0.50 

14 

3.0 

3.17 

5.75 

18 

3.0 

2.96 

1.50 

6 

5.0 

4.87 

2.65 

10 

5.0 

5.04 

0.80 

14 

5.0 

4.81 

3.85 

18 

5.0 

5.07 

1.40 

6 

7.0 

6.80 

2.82 

10 

7.0 

7.38 

5.43 

14 

7.0 

7.10 

1.46 

18 

7.0 

7.42 

5.96 
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Wave  direction 


Two  wave  directions  were  chosen  to  reproduce  maximum  vessel  excursions. 
Vessel  pitch  and  heave  motion  (see  Figure  11)  are  most  strongly  excited  when 
waves  encounter  the  vessel  from  either  front  or  back.  Waves  approaching 
directly  along  the  channel  would  reproduce  this  condition.  When  waves 
approach  broadside  or  at  90  deg  from  the  long  axis  of  the  vessel,  roll  motion  is 
most  strongly  excited.  This  would  imply  that  waves  are  traveling  perpendicular 
to  the  beach.  Refraction  makes  this  wave  approach  unlikely.  However,  a 
combined  roll/pitch  motion  is  strongly  excited  when  waves  travel  approximately 
20  deg  from  the  long  axis  of  the  vessel.  According  to  harbor  pilots,  waves  from 
the  southwest  generate  considerable  vessel  roll.  The  wavemaker  was  positioned 
so  that  waves  are  generated  20  deg  south  of  the  entrance  channel  orientation. 

The  plunger  wavemaker  is  a  unidirectional  machine;  thus,  to  generate  both  wave 
directions,  the  wavemaker  had  to  be  physically  moved.  The  two  wavemaker 
positions  are  shown  in  Figures  18, 19,  and  20. 


Modifications  to  Existing  Transition 

Transition  location 

Presently  the  Barbers  Point  entrance  channel  is  dredged  to  -12.8  m  (-42  ft) 
mllw  and  the  harbor  to  -11.6  m  (-38  ft)  mllw.  The  transition  between  these  two 
depths  is  located  at  the  intersection  of  the  shoreline  with  the  entrance  channel. 
During  early  stages  of  this  study,  it  became  apparent  that  the  transition  would  be 
a  point  of  possible  groundings.  The  waves  are  fairly  energetic  in  the  channel  at 
the  shoreline.  It  was  decided  to  move  the  transition  from  the  shoreline  location 
to  where  the  channel  opens  up  to  the  harbor  (Figure  24).  This  decision  was 
based  upon  two  reasons:  decreased  wave  energy  at  the  second  transition  and 
vessel  shear  (described  next). 

When  a  vessel  transits  over  a  bottom  that  goes  from  deeper  to  shallower 
water,  the  vessel  will  shear  if  the  water  is  shallow  enough.  This  shear  is  caused 
by  increased  boundary  layer  friction  between  the  vessel  moving  through  the 
water  and  the  seafloor  bottom.  The  portion  of  the  boat  in  deeper  water  (stem) 
feels  less  friction  and  wants  to  travel  more  quickly  than  that  part  of  the  vessel  in 
shallower  water  (bow).  The  vessel’s  stem  will  try  to  “fishtail”  or  “jackknife.” 

By  moving  the  transition  toward  the  opening  of  the  harbor,  the  pilot  has  more 
room  to  correct  the  vessel  position. 

The  transition  was  moved  to  the  new  location,  and  all  data  were  collected 
with  the  transition  at  this  new  location.  The  existing  difference  in  the  depth 
between  the  entrance  channel  and  the  harbor  of  1.2  m  (4  ft)  was  used  for  the  first 
series  of  experiments. 
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Figure  24.  Location  of  initial  and  proposed  transition  between  entrance  channel 
and  harbor 


Difference  in  entrance  channel  and  harbor  depths 

Although  there  is  only  0.6-m  (2-ft)  (  State  HDOT  criteria)  or  1.2-m  (4-ft) 
(Corps  POD  criteria)  underkeel  clearance  in  the  harbor,  there  is  1.8-m  (6-ft) 
(State  HDOT  criteria)  or  2.4-m  (8-ft)  (Corps  POD  criteria)  underkeel  clearance 
in  the  entrance  channel.  The  entrance  channel  is  approximately  1,070  m 
(3,500  ft)  long  and  140  m  (450  ft)  wide.  There  is  a  potential  for  significant 
savings  if  depth  in  the  entrance  channel  can  be  minimized.  During  the  first 
series  of  experiments,  it  was  determined  that  the  entrance  channel  could  be  made 
shallower.  Thus,  navigation  experiments  were  also  conducted  with  a  shallower 
entrance-channel  depth  and  a  smaller  depth  difference  at  the  transition.  A  0.6-m 
(2-ft)  channel  transition  was  also  modeled.  By  adjusting  the  water  depths, 
underkeel  clearance  values  for  the  entrance  channel  and  harbor,  respectively,  of 
1.8  m  (6  ft)  and  1.2  m  (4  ft),  1.5  m  (5  ft)  and  0.9  m  (3  ft),  and  1.2  m  (  4  ft)  and 
0.6  m  (2  ft)  were  evaluated. 


Design  Vessel  Draft 

and  Water-Level  Combinations 

As  previously  discussed,  there  is  a  difference  between  the  HDOT  and  the 
Corps  POD  underkeel  clearance  criteria  in  the  harbor.  HDOT  uses  a  0.6-m  (2-ft) 
clearance,  while  POD  uses  a  1.2-m  (4-ft)  underkeel  clearance  value.  It  was 
decided  that  the  design  vessels  would  be  drafted  to  reproduce  both  the  HDOT 
and  the  POD  underkeel  clearance  criteria  for  the  selected  water  depths. 
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A  total  of  1,280  inbound  and  outbound  runs  were  conducted  with  the  two 
design  vessels.  Table  6  summarizes  the  vessel  draft/water  depth  combinations 
evaluated  by  this  physical  model  study. 


Table  6 

Vessel  Draft  and  Water  Depth  Combinations  Evaluated  in  the  Physical  Model 

Vessel 

Water  Depth,  ft 

Underkeel  Clearance,  ft 

Number  of  Runs  | 

Description 

Draft,  ft 

Channel 

Harbor 

Channel 

Harbor 

Incident  Wave  Directions 

Total 

0 

20 

APL  C9-class 
containership 

35 

41 

37 

6 

2 

48 

48 

96 

35 

43 

39 

8 

4 

48 

48 

96 

Modified  Bunga 
Saga  Empat  bulk- 
cargo  carrier 

38 

44 

40 

6 

2 

48 

48 

96 

36 

44 

40 

8 

4 

48 

48 

96 

36 

42 

38 

6 

2 

48 

48 

96 

34 

42 

38 

8 

4 

48 

48 

96 

39 

47 

43 

8 

4 

48 

48 

96 

41 

47 

43 

6 

2 

48 

48 

96 

41 

49 

45 

8 

4 

48 

48 

96 

43 

49 

45 

6 

2 

48 

48 

96 

41 

47 

45 

6 

4 

16 

16 

32 

41 

46 

44 

5 

3 

16 

16 

32 

41 

45 

43 

4 

2 

16 

16 

32 

43 

49 

47 

6 

4 

16 

16 

32 

43 

48 

46 

5 

3 

16 

16 

32 

43 

47 

45 

4 

2 

16 

16 

32 

39 

45 

43 

6 

4 

16 

16 

32 

39 

44 

42 

5 

3 

16 

16 

32 

39 

43 

41 

4 

2 

16 

16 

32 

42 

47 

45 

5 

3 

16 

16 

32 

Total 

1,280 

Vessel  Speed 

Harbor  pilots  provided  typical  vessel  transit  speeds  into  and  from  the  harbor. 
They  noted  that  vessel  speeds  on  the  inbound  transit  varied  from  7  knots  at  the 
channel  entrance  to  2  knots  in  the  harbor.  On  the  outbound  transit,  vessel  speed 
increased  from  2  knots  in  the  harbor  up  to  10  knots  at  the  channel  entrance 
during  rough  conditions  and  strong  currents.  Vessel  speeds  of  approximately 
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iments  around  the  transition  as  well  as  in  the 


5  Physical  Model  Results 
and  Conclusions 


Over  nine  gigabytes  of  data  were  collected  during  the  approximately 
1,000  transits  into  and  out  of  the  harbor  for  the  six  virtual  markers  on  the 
modified  Bunga  Saga  Empat  bulk-cargo  carrier  and  the  four  virtual  markers  on 
the  APL  C9-class  containership.  All  the  data  are  not  presented  here,  but 
summarized  data  are  provided  in  tabular  and  graphical  form. 

Early  in  the  study,  it  became  apparent  that  the  area  around  the  transition 
could  be  a  region  of  maximum  vertical  excursion  of  the  vessel.  The  majority  of 
the  data-collection  efforts  were  focused  on  this  area.  Accurately  locating  the 
channel  depth  was  critically  important  since  the  channel  is  945  m  (3,100  ft)  long. 
A  concentrated  effort  with  large  waves  was  also  conducted  in  the  entrance 
channel. 


Data  Analysis  Techniques 

During  experimental  runs,  data  were  collected  at  60  Hz.  At  the  end  of  the 
data-collection  effort,  video  images  of  the  reflecting  balls  were  analyzed,  and  the 
x-,  y-,  z-  coordinates  of  the  six  balls  located  on  the  model  vessel  and  the  balls 
located  on  the  floor  around  the  channel  were  identified  manually.  After  each  of 
the  reflecting  balls  were  identified,  the  computer  calculated  the  x-,  y-,  z- 
coordinate  time-history  of  each  reflecting  ball  using  the  arbitrarily  defined 
coordinate  system.  Then,  the  x-,  y-,  z-  coordinates  were  transposed  into  the 
model  coordinate  system  using  the  reflecting  balls  located  around  the  channel. 
Assuming  the  model  vessel  was  rigid,  the  locations  of  the  virtual  markers  were 
calculated  at  each  time-step.  The  vessel  motion  then  was  compared  with  a  static 
snapshot  of  the  model  vessel  in  the  center  of  the  sampling  area.  After  these 
calculations  were  made,  the  underkeel  clearance  of  the  moving  vessel  could  be 
determined  since  the  static  underkeel  clearance  was  measured  with  a  caliper. 

Once  the  locations  of  the  virtual  markers  were  known,  the  underkeel 
clearance  was  calculated,  and  the  data  were  sorted  into  0.3-m  (1-ft)  bins 
representing  the  keel  distance  above  the  bottom  for  each  of  the  virtual  markers 
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for  all  the  navigation  runs.  Grouping  the  data  into  depth  bins  instead  of 
presenting  time-histories  allowed  the  data  to  be  further  summarized  and 
presented  in  tabular  form.  Instead  of  having  six  different  time-histories  of  the 
location  of  the  model  vessel  keel,  the  data  for  all  the  virtual  markers  were 
grouped  together.  This  decreased  the  number  of  separate  channels  of  data  from 
six  to  one  for  the  modified  Bunga  Saga  Empat  and  from  four  to  one  for  the  APL 
C9-class  vessel  for  each  experimental  condition.  Also,  the  data  for  the  inbound 
and  outbound  runs  were  combined  into  one  data  set  for  each  wave  condition. 

For  simplicity,  the  percent  occurrence  in  each  bin  is  presented  instead  of 
reporting  the  number  of  occurrences  of  each  virtual  marker  in  each  depth  bin. 
The  data  have  been  scaled  into  prototype  distances,  with  the  results  for  the 
modified  Bunga  Saga  Empat  presented  in  Appendixes  A  and  B  for  Corps  POD 
and  State  HDOT  criteria,  respectively.  The  data  are  presented  by  the  vessel 
draft/entrance  channel  depth/harbor  depth  combinations.  Appendixes  C  and  D 
show  the  results  for  the  APL  C9-class  vessel  for  Corps  POD  and  State  HDOT 
criteria,  respectively. 

Figure  25  shows  the  time-history  of  the  virtual  marker  trajectories  for  one 
wave  case  for  the  three  virtual  markers  at  the  front  of  the  modified  Bunga  Saga 
Empat  model  vessel  (location  of  virtual  markers  shown  in  Figures  22  and  23). 
The  left  virtual  marker  is  defined  as  the  port  side,  and  the  right  virtual  marker  is 
the  starboard  side  of  the  model  vessel.  Figure  25  shows  the  results  from  the 
second  camera  arrangement  with  the  sampling  area  entirely  in  the  channel.  The 
data  are  unedited,  and  some  spikes  in  the  data  are  visible.  These  would  be 
eliminated  during  data  quality  check.  The  x-axis  in  Figure  25  shows  the  location 
along  the  channel,  and  the  y-axis  indicates  the  distance  from  the  bottom. 


Results  for  the  Modified  Bunga 
Saga  Empat  Bulk-Cargo  Carrier 

Tabulated  results  are  organized  by  the  underkeel  clearance  criteria  used 
during  the  experiment  and  by  the  transition  difference.  Table  7  shows  percent 
occurrence  of  keel  motion  for  the  Corps  POD  underkeel  clearance  criteria,  while 
Appendix  A  shows  the  actual  number  of  occurrences  of  keel  motions  for  the  first 
camera  arrangement,  the  proposed  transition  location,  and  a  1.2-m  (4-ft) 
transition  difference.  Table  8  and  Appendix  B  show  corresponding  results  for 
the  State  HDOT  underkeel  clearance  criteria.  The  depth  bins  are  organized  in 
0.3-m  (1-ft)  increments  where  “>8”  indicates  that  the  vessel  keel  is  2.4  m  (8  ft) 
above  the  bottom.  The  bins  below  “>8”  indicate  that  the  keel  is  identified  by  the 
virtual  marker  trajectories  passed  in  this  zone.  For  example,  for  the  36/44/0  case 
(vessel  draft  of  36  ft,  water  depth  in  the  channel  of  44  ft,  and  waves  arriving 
directly  along  the  entrance  channel)  for  the  “>2”  ft  depth  bin,  the  keel  was 
between  0.9  m  (3  ft)  and  0.6  m  (2  ft)  from  the  bottom  0.367  percent  of  the  time 
in  the  harbor.  The  results  for  the  harbor  are  shown  above  the  results  for  the 
channel.  The  last  column  indicates  the  average  number  of  occurrences  of  the 
keel  in  each  depth  bin  over  the  numerous  vessel  draft/water  depth  and  wave 
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Virutal  marker  elevation  (model  in)  Vi rutal  marker  elevation  (model  in)  Vi rutal  marker  elevation  (model  in) 


Figure  25.  Unedited  track  of  three  virtual  markers  for  the  modified  Bunga  Saga 
Empat  bulk-cargo  carrier  in  the  channel,  1 .5-m  (5-ft),  14-sec 
significant  wave 
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direction  combinations.  The  average  column  indicates  that  the  Corps  POD 
ntena  provide  sufficient  underkeel  clearance  in  the  channel  and  harbor  For 
example,  the  vessel  is  between  0  6  m  (2  frt  and  tn  m  n  m  f  ,,  '  or 

0.001  nercenf  of  tho  t,w  ™  K  }  d  °’3  m  (1  from  the  bo«om 

near  the I™  1?  The  aveIaSe  ™*»]  speed  in  knots  is  also  provided 

ear  me  top  of  the  table  for  each  condition  studied. 

jSSriS-aKSis-' 

clearance  distance  is  not  adequate,  and  the  1.2-m  (4-ft)  clearance  distance  is 

overly  conservative.  Both  the  1.8-m  (6-ft)  and  2  4  m  fR  m  , 

the  channel  appear  to  be  conservative  }  (M°  ClearaDce  d,slances  in 

0  6  ftTal  fT  conducted  with  a  modified  transition  depth  of 

?*o  ” <? “  no  « i'2 m <4 ft)'  Underkeel cleara”“s of  1.2 m (4 fO  and 

tobor^h  m,(  #)  “d  1S  m  (5  *>•  aad  0  6  ”  <2  *)  »"d  1-2  m  (4  ft)  in  the 
harbor  and  channel,  respectively,  were  evaluated.  The  percent  occurences  of 

iTbleT T 1 2^  SftT  °fd1  have  bee”  rompiled  and  Me  shown  in 

“  «-t  r*.  SSSST5S 

.  .  ^  .  )  an<^  1-2-m  (4-  ft)  underkeel  clearances  in  the  harbor  and 

hannel,  respectively.  Appendix  E  shows  the  actual  number  of  keel  motion*  f 
the  appropriate  corresponding  conditions.  (The  43-ft  vessel  draft  47  ft  1m  * 
channe,  depth,  45-ft  harbor  depth,  waves  0-deg 

ft  transition  at  proposed  location  data  were  inadvertently  omitted  from  this  ’ 
appendix.)  Results  indicate  that  the  1.2-m  (4-ft)  and  1  8  m  TfTt  n  , 
channel  underkeel  clearances,  respectively,' ^  ^ 

Data  were  also  collected  in  the  entrance  channel  with  the  second  camera 
arrangement.  During  this  data  collection,  wave  heights  were  run  multmle  t,W* 

achT/rhT  T**  Valid‘,y'  The 

each,  and  the  results  were  summarized  for  each  different  replication  The 
percent  occurrence  of  keel  motion  for  these  conditions  are  shown  in  Table  12 
The  modified  Bunga  Saga  Empat  was  drafted  so  that  there  Z  a 1 5-m  5  ft) 
channel  underkeel  clearance.  The  depth  in  the  channels  was  14  3  m  (47ft  L 
the  incident  wave  direction  was  zero.  Appendix  F  shows  the  actual  nnmfJr  nf 
occurrences  of  keel  motions  for  the  appropriate  corresponding  conditions 
Results  indicate  that  a  1.5-m  (5-ft)  underkeel  clearance"*  the^hannd  should 
provide  adequate  clearance  for  the  design  vessel. 

Results  for  the  APL  C9-Class  Containership 

BunraeSa»aPEmn!',e  C9'ClaSS  VeSsel  is  more  s,rMm]i“ed  than  the  modified 
ga  Saga  Empat,  and  experiments  were  conducted  with  this  model  vessel  to 
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Tab,e  9  *„»!  Motion  Through  Different  Depth  Bins  for  Modified  Bungs 

Percent  Occurrence  of  Keel  Motion  Throug  c|earance  criteria  for  Varying 

Saga  Empat  with  (Underkeel 

--  are  1.8  m  (6  «,  and  1.2  ai,a  «), 

41/47/0  U3,4S„  39/45/0  \«~  1^ 

channel/wave  direction - - - - -  .  R 

- - -  44  45  4.3  4.2  4.2  j  4.6  J - 

Speed,  knots _ L_ - L— - '  _ 


Depth  bins,  ft 

>8.0 _ 

>7.0 _ 

>6.0 

>5.0 

>4.0 

>3.0 

>2.0 

>1.0 _ 

>0.0 

Bottom 

Total 

Depth  bins,  ft 
>8.0 


Percent  Occurrence  in  Each  Depth  Bin  over  the  Harbor 


53.086  36.198 

44.993  62.317 


0.000 

0.000 

0.000 

0.000 

Average 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0  065 

0.000 

0.153 

0.000 

0.339 

0.000 

0.030 

0.000 

0.177 

7.066 

84.243 

28.180 

69.742 

39.075 

59.639 

7.938 

82.358 

25.473 

70.192 

8.626 

1.925 

0.946 

9.578 

4.131 

0.000 

0.000 

0.000 

0.096 

0.027 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Percent  Occurrence  in  Each  Depth  Bin  over  the  Entrance  Channel_ 


42.416  34.331 

56.180  64.499 


0.099 

6.130 

86.083 

7.668 


0.020 


0.000 


0.000 


0.000 


0.000 


0.307 

28.605 

69.535 

1.433 


0.057 


0.000 


0.000 


0.000 


0.000 


0.584 
37.976 
59.584 
'  1.844 


0.012 


0.000 


0.000 


0.000 


0.000 


0.097 

7.227 

85.013 

7.644 


0.019 


0.000 


0.000 


0.000 


0.000 


Average 

0.321 

23.426 

72.238 

3.974 


0.026 


0.000 


0.000 


0.000 


0.000 


Navigation^ experiments  with^ toe^^S-class  vessel  were  conducted  with  the 

J2  “subs  for  the  four  coupons  studied,  prevously  shown  in  Table  6. 
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Table  10 

Percent  Occurrence  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified  Bunga 
Saga  Empat  with  Modified  Transition  for  Underkeel  Clearance  Criteria  for  Varying 
Vessel  Draft/Water  Depth  in  Harbor/Wave  Direction  Occurrences  (Underkeel  clearances 
in  the  entrance  channel  and  harbor  are  1.5  m  (5  ft)  and  0.9  m  (3  ft),  respectively) 


Vessel  draft/water  depth  in  channel/wave  direction 


peed,  knots 


Depth  bins,  ft 


>8.0 


>7.0 


>6.0 


>5.0 


>4.0 


>3.0 


>1.0 


>0.0 


Bottom 


otal 


Depth  bins,  ft 


>8.0 


>7.0 


>6.0 


>5.0 


>4.0 


>3.0 


-2.0 


>1.0 


>0.0 


41/46/0  43/48/0  39/44/0  41  /46/20 


.5 


.2 

4 

.3 

4 

.3 

4 

43/48/20 


4.1 


39/44/20 


Percent  Occurrence  in  Each  Depth  Bin  over  the  Harbor 


Average 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.201 

0.000 

0.014 

0.222 

0.366 

0.055 

0.196 

29.116 

35.075 

3.019 

27.655 

40.219 

5.084 

22.792 

66.064 

62.687 

92.248 

69.397 

57.636 

88.135 

73.344 

4.618 

2.239 

4.720 

2.703 

1.777 

6.717 

3.658 

0.000 

0.000 

0.000 

0.023 

0.003 

0.008 

0.010 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Percent  Occurrence  in  Each  Depth  Bin  over  the  Entrance  Channel 


0.000  0.000 


1.515  0.231 


33.333  23.095 


0.000  0.016 


0.000  0.254 


3.651  28.557 


62.626  76.212  91.074  69.259 


2.525  0.462  5.269  1.910 


0.000  0.000  0.006  0.004 


0.000  0.000 


0.000  0.000 


0.000  0.000  0.000  0.000 


0.000  0.000  0.000  0.000 


Average 


0.000  0.000  0.004 


0.291  0.074  0.166 


39.226  5.982  19.941 


58.785  88.022  76.238 


1.650  5.852  3.617 


0.048  0.067  0.033 


0.000  0.004  0.001 


0.000  0.000  0.000 


0.000  0.000  0.000 


Table  13  shows,  for  both  the  Corps  POD  and  State  HDOT  underkeel  clearance 
criteria,  the  percent  occurrence  of  keel  motion,  while  Appendixes  C  and  D  show 
the  actual  number  of  occurrences  of  keel  motions  for  these  two  underkeel 
clearance  criteria  for  appropriate  corresponding  conditions.  The  percent 
occurrences  of  all  the  virtual  markers  for  all  the  wave  conditions  have  been 
included.  These  results  indicate  that  using  the  State  HDOT  criteria  could  likely 
result  in  groundings. 
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Vessel  draft/water  depth  in  channel/wave  direction 

41/45/0 

43/47/0 

39/43/0 

Speed,  knots 

3.3 

3.8 

3 

Table  11 

Percent  Occurrence  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified  Bunga 
aga  Empat  with  Modified  Transition  for  Underkeel  Clearance  Criteria  for  Varying 
Vessel  Draft/Water  Depth  in  Channel/Wave  Direction  Occurrences  (Underkeel 
clearances  in  the  entrance  channel  and  harbor  are  1.2  m  (4  ft)  and  0.6  m  (2  ft), 


43/47/20 

39/43/20 

3.5 

2.6 

Percent  Occurrence  in  Each  Depth  Bin  over  the  Harbor 


Depth  bins,  ft 


Average 


>8.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

>7.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

>6.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

>5.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

>4.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

>3.0 

0.271 

0.000 

0.490 

0.295 

0.622 

0.051 

0.466 

>2.0 

28.424 

66.667 

2.941 

28.659 

40.590 

6.339 

6.206 

>1.0 

67.894 

33.333 

85.907 

69.287 

57.101 

88.609 

83.800 

>0.0 

3.411 

0.000 

10.662 

1.753 

1.687 

4.954 

9.525 

Bottom 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

Percent  Occurrence  in  Each  Depth  Bin  over  the  Entrance  Channel 


Depth  bins,  ft 


>8.0 


>7.0 


>6.0 


>5.0 


Average 


0.000 

0.000 

0.000 

0.012 

0.000 

0.000 

0.003 

0.000 

0.000 

0.000 

0.024 

0.005 

0.000 

0.006 

0.000 

0.000 

0.000 

0.131 

0.372 

0.019 

0.102 

28.571 

66.667 

0.000 

28.715 

39.297 

1.695 

14.189 

69.870 

33.333 

97.333 

70.068 

59.107 

95.870 

83.717 

1.558 

0.000 

2.667 

1.050 

1.219 

2.415 

1.983 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

42 


Chapter  5  Physical  Model  Results  and  Conclusions 


Table  12 

Percent  Occurrence  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified  Bunga 
Saga  Empat  in  Entrance  Channel  for  1.5-m  (5-ft)  Underkeel  Clearance  (The  four  largest 
waves  were  run  eight  times  each) 


[Depth  bins,  ft 


>8.0 


>7.0 


>6.0 


>5.0 


>4.0 


b>3.0 


>2.0 


fc>1.0 


>0.0 


Bottom 


Runl 


0.009 


0.009 


0.678 


48.965 


49.004 


1.331 


0.004 


0.000 


0.000 


0.000 


Run2 


0.000 


0.013 


0.471 


49.192 


48.277 


2.027 


0.020 


0.000 


0.000 


0.000 


Run3 


0.000 


0.013 


1.033 


48.079 


49.454 


1.309 


0.112 


0.000 


0.000 


0.000 


Run4 


0.006 


0.012 


0.744 


47.045 


51 .729 


0.463 


0.000 


0.000 


0.000 


0.000 


Run5 


0.000 


0.000 


0.650 


46.457 


52.286 


0.599 


0.007 


0.000 


0.000 


0.000 


Run6 


0.000 


0.008 


0.721 


49.598 


49.109 


0.564 


0.000 


0.000 


0.000 


0.000 


Run7 


0.000 


0.000 


0.801 


48.537 


49.148 


1.455 


0.017 


0.043 


0.000 


0.000 


Run8 


0.050 


0.010 


0.915 


47.141 


50.920 


0.955 


0.010 


0.000 


0.000 


0.000 


Average 


0.008 


0.008 


0.752 


48.127 


49.991 


1.088 


0.021 


0.005 


0.000 


0.000 


Study  Conclusions 


Based  on  analysis  of  the  physical  model  data,  the  following  conclusions  and 
recommendations  are  deduced: 


a.  The  State  HDOT  underkeel  criterion  in  the  harbor  does  not  allow  enough 
clearance  for  the  design  vessels  to  transit  into  the  harbor  without  possible 
groundings. 

b.  The  State  HDOT  and  Corps  POD  criteria  in  the  channel  are  both 
conservative,  and  no  groundings  or  near  groundings  were  found. 

c.  The  recommended  underkeel  clearance  is  0.9  m  (3  ft)in  the  harbor  and 
1.5  m  (5  ft)  in  the  channel  for  the  design  vessels  studied  and  for  waves 
less  than  2.1  m  (7  ft)  in  height. 

d.  The  transition  should  be  moved  to  the  opening  of  the  harbor  basin.  At 
the  proposed  transition  location,  there  is  less  wave  energy;  if  vessel  shear 
occurs,  the  harbor  pilot  would  have  more  room  to  react/correct.  Model 
navigation  study  data  support  a  0.6-m  (2-ft)  transition. 
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Table  13 

Percent  Occurrence  of  Keel  Motion  Through  Different  Depth  Bins  for  APL 
C9-Class  Containership  for  Corps  POD  and  State  HDOT  Underkeel  Clearance 
Criteria  with  Proposed  Transition  Location  and  4-ft  Depth  Difference  at  the 
ransition 


Corps  POD 

State  HDOT 

35/43/0  deg 

35/43/20  deg 

35/41/0  deg 

35/41/20  deg 

Average 

Average 

Average 

Averages 

peed,  knots 


Depth  bins,  ft 


>8.0 


>2.0 


>1.0 


>0.0 


Bottom 


Depth  bins,  ft 


>2.0 


>1.0 


>0.0 


Bottom 


Percent  Occurrence  in  Each  Depth  Bin  over  the  Harbor 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.036 

0.002 

0.000 

0.000 
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Appendix  A 

U.S.  Army  Engineer  Division, 
Pacific  Ocean  (POD), 

Underkeel  Clearance  Criteria 
Results  for  the  Modified  Bunga 
Saga  Empat  Bulk-Cargo  Carrier 


The  following  abbreviation  codes  have  been  used  throughout  Appendixes  A 
through  F  to  describe  the  wave  conditions,  vessel  direction,  and  run  number.  In 
general,  two  runs  were  made  for  each  condition. 


Abbreviation  Code 

Wave  Condition 

Vessel  Direction 

Run  Number 

lil 

1 

Into  harbor 

1 

1i2 

1 

Into  harbor 

2 

1  ul 

1 

Out  of  harbor 

1 

1u2 

1 

Out  of  harbor 

2 

2i1 

2 

Into  harbor 

1 

2i2 

2 

Into  harbor 

2 

2u1 

2 

Out  of  harbor 

1 

2u2 

2 

Out  of  harbor 

2 

3i1 

3 

Into  harbor 

1 

3i2 

3 

Into  harbor 

2 

3u1 

3 

Out  of  harbor 

1 

3u2 

3 

Out  of  harbor 

2 

4i1 

4 

Into  harbor 

1 

4i2 

4 

Into  harbor 

2 

4u1 

4 

Out  of  harbor 

1 

(continued) 
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The  average  vessel  speed  is  given  as  a  vector  quantity  based  upon  the  coordinate 
convention  shown  in  Figures  19  and  20  of  the  main  text.  For  the  average  speed 
calculated  for  a  given  condition,  the  absolute  value  of  the  speed  was  used.  The 
number  of  occurrences  is  calculated  by  summing  each  occurrence  of  each  virtual 
marker  at  each  instant  in  time.  This  is  different  from  the  body  of  the  report 
where  the  percent  of  occurrence  is  reported. 
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Saga  Empat  Bulk-Cargo  Carrier 


See  Appendix  A  for  description  of  tables. 
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Location 
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Table  El  2 

Number  of  Occurrences  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified 
Bunga  Saga  Empat,  39-ft  Vessel  Draft,  45-ft  Entrance-Channel  Depth,  43-ft  Harbor 
Depth,  Waves  O-deg  Orientation  to  Entrance  Channel,  2-ft  Transition  at  Proposed 
Location 
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Table  El  3 

Number  of  Occurrences  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified 
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Location 
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Table  El  4 

Number  of  Occurrences  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified 
Bunga  Saga  Empat,  39-ft  Vessel  Draft,  a  43-ft  Entrance-Channel  Depth,  41 -ft  Harbor 
Depth,  Waves  0-deg  Orientation  to  Entrance  Channel,  2-ft  Transition  at  Proposed 
Location 
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Number  of  Occurrences  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified 
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Location 
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Appendix  F 

Results  for  the  Modified  Bunga 
Saga  Empat  Bulk-Cargo  Carrier 
in  the  Entrance  Channel 


See  Appendix  A  for  description  of  tables. 
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Table  F5 

Number  of  Occurrences  of  Keel  Motion  Through  Different  Depth  Bins  for  Modified  Bunga  Saga  Empat,  42-ft  Vessel 
Draft,  47-ft  Entrance-Channel  Depth,  Waves  O-deg  Orientation  to  Entrance  Channel,  Run  5  in  Entrance  Channel _ 

Code  9(5  9u5  ai5  au5  bi5  bu5  ci5  cu5  Total 
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Appendix  G 

Plan  4c  Impact  on  Barge  Basin 
Wave  Climate1 


Harbor  Description 

The  harbor  complex  presently  consists  of  an  entrance  channel,  deep-draft 
harbor,  barge  basin,  and  a  resort  marina.  The  entrance  channel  is  140  m  (450  ft) 
wide,  945  m  (3,100  ft)  long,  and  12.8  m  (42  ft)  deep.  The  deep-draft  harbor 
basin  is  11.6  m  (38  ft)  deep.  The  barge  basin,  located  just  seaward  of  the  harbor 
on  the  south  side  of  the  entrance  channel,  is  relatively  exposed  to  incident  wave 
energy.  It  is  67  m  (220  ft)  wide,  396  m  (1,300  ft)  long,  and  7.0  m  (23  ft)  deep. 
The  entrance  channel  is  aligned  approximately  with  the  225-deg  azimuth  (from 
SW  to  NE  orientation).  The  existing  harbor  configuration  is  designated  as 
Plan  la,  as  in  previous  studies  (e.g.,  Briggs  et  al.  19942). 

Planned  modifications  to  Barbers  Point  Harbor  include  a  deeper  entrance 
channel  and  a  rubble-mound  jetty  on  the  north  side  of  the  entrance  (Plan  4c). 

The  modifications  will  impact  wave  and  swell  conditions  in  the  relatively 
exposed  barge  basin.  The  Pacific  Ocean  Division  requested  the  U.S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  to  assess  the  impact  of  Plan  4c 
on  barge  basin  wave  climate.  WES  used  wave  climate  data  from  the  offshore  Sxy 
array  gauge  and  the  numerical  model  HARBD  to  assess  overall  barge  basin  wave 
climate  in  both  the  existing  and  Plan  4c  harbors.  Results  from  previous  physical 
model  storm  wave  simulations  were  also  reviewed  to  assess  barge  basin  response 
to  storms. 

The  proposed  Plan  4c  includes  the  following:  (a)  flare  of  the  outer  305  m 
(1,000  ft)  of  the  entrance  channel  from  a  width  of  140  m  (450  ft)  to  229  m 
(750  ft),  (b)  140-m-  (450-ft-)  long  shore-connected  rubble-mound  jetty  on  the 


1  Appendix  G  was  prepared  by  Dr.  Edward  F.  Thompson  and  Messrs.  Michael  J.  Briggs  and 
Doyle  L.  Jones. 

2  Briggs,  M.,  Lillycrop,  L.,  Harkins,  G.,  Thompson,  E.,  and  Green,  D.  (1994).  “Physical  and 
numerical  model  studies  of  Barbers  Point  Harbor,”  Technical  Report  CERC-94-14,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 
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G1 


north  side  of  the  entrance  channel,  (c)  channel  and  harbor  deepened  to  14.9  m 
(49  ft)  and  13.7  m  (45  ft),  respectively,  and  (d)  335-m  (1,100-ft)  by  335-m 
(1,100-ft)  harbor  expansion. 


Incident  Wave  Climate 

Barbers  Point  Harbor  is  subject  to  waves  approaching  the  Hawaiian  Islands 
from  the  northwest  and  the  southwest.  Deepwater  wave  approach  from  the  west 
occurs  rarely  during  Kona  (local)  storms.  The  island  of  Oahu  shelters  the  harbor 
from  waves  generated  by  easterly  trade  winds.  The  largest  waves,  occurring 
during  winter  months,  are  swell  because  of  storms  in  the  northwest  Pacific. 

Most  waves  measured  at  the  SX7  gauge  location  outside  the  harbor  complex  have 
peak  periods  between  9  and  15  sec  with  significant  wave  heights  ranging  from 
0.3  m  (1  ft)  to  0.9  m  (3  ft). 

Wave  parameters  for  analysis  from  the  Sxy  gauge  data  were  obtained  from  the 
Scripps  Institution  of  Oceanography.  Data  covered  the  period  September  1986 
to  October  1989.  Parameters  were  used  to  develop  percent  occurrence  tables  of 
significant  wave  height,  peak  wave  period,  and  dominant  wave  direction.  Wave 
height  intervals  are  0.2  m,  and  wave-direction  intervals  are  10  deg. 


Wave  Climate  in  Barge  Basin 

The  numerical  model  HARBD  was  used  to  calculate  wave  transformation 
between  incident  waves  (Sxy  gauge  location)  and  waves  in  the  barge  basin  for 
both  the  existing  harbor  and  Plan  4c.  Model  grids  used  in  previous  long  wave 
studies  (e.g.,  Briggs  et  al.  1994 x)  were  adapted  for  this  wave  and  swell 
application.  Wave  periods  from  9  to  20  sec,  in  1-sec  intervals,  and  wave 
directions  from  190-  to  280-deg  azimuth,  in  10-deg  intervals,  were  run  in  the 
model. 

HARBD  results  were  further  processed  to  provide  directional  spectral 
estimates  of  wave  transformation  at  selected  points  around  the  harbor.  The 
procedure  is  described  by  Thompson  et  al.  (1996).1 2  Transformation  cases  were 
matched  to  wave  period  and  direction  intervals  used  for  the  wave  gauge  climate 
summaries. 

Wave  climate  percent  occurrence  tables  from  the  gauge  were  then 
transformed  into  percent  occurrence  tables  at  two  points  along  the  approach  to 


1  Briggs,  M.,  Lillycrop,  L.,  Harkins,  G.,  Thompson,  E.,  and  Green,  D.  (1994).  “Physical  and 
numerical  model  studies  of  Barbers  Point  Harbor,”  Technical  Report  CERC-94-14,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

2  Thompson,  E.,  Hadley,  L.,  Brandon,  W.,  McGehee,  D.,  and  Hubertz,  J.  (1996).  “Wave 
response  of  Kahului  Harbor,  Maui,  Hawaii,”  Technical  Report  CERC-96-11,  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 
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the  harbor  (channel  center  about  500  ft  seaward  of  the  jetty  tip  and  channel 
center  at  the  beginning  of  the  confined  channel  in  line  with  the  landward  jetty 
end)  and  three  points  in  the  barge  basin  (north  part  of  the  barge  basin  and  the 
east  and  south  barge  basin  corners;  Figure  Gl).  Numerical  model  output  station 
numbers  are  also  shown  in  Figure  Gl. 


Figure  Gl .  Numerical  model  output  station  locations 

Mean  significant  wave  height,  Hs,  and  maximum  significant  wave  height, 
Hs,max  ,  were  also  calculated  for  each  point  (Figures  G2  throughG7,  and 
Table  Gl).  Barge  basin  Hs  variation  with  direction  shows  a  maximum  for  waves 
approaching  the  entrance  from  the  west  (270  deg).  Hs  for  waves  from  the  west 
in  Plan  4c  is  slightly  lower  than  in  the  existing  harbor,  indicating  that  the  jetty  is 
providing  some  shelter  to  the  barge  basin.  For  waves  coming  from  south  of 
approximately  240  deg,  Hs  in  Plan  4c  exceeds  H s  in  Plan  la,  indicating  the 
deeper  channel  impacts  barge  basin  wave  conditions.  Overall  barge  basin  Hs  is 
the  same  for  both  Plan  4c  and  existing  harbors. 

The  Hs  max  is  more  variable,  but  tends  to  follow  a  similar  trend.  Maximum 
values  in  the  barge  basin  comers  are  lower  in  Plan  4c  than  in  Plan  la  since  the 
maximum  wave  approached  from  the  west  was  partially  blocked  by  the  jetty. 

The  standard  Corps  criterion  for  an  acceptable  mooring  area  is  that  Hs  should 
not  exceed  0.3  m  (1  ft)  more  than  10  percent  of  the  time.  The  barge  basin  was 
evaluated  relative  to  this  criterion  (Table  G2).  Comer  areas  met  the  criterion  in 


Appendix  G  Plan  4c  Impact  on  Barge  Basin  Wave  Climate 


G3 


Mean  Wave  Height  (ft) 
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Figure  G2.  Mean  Hs  from  HARBD,  barge  basin  south  corner,  Station  7 
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Figure  G3.  Mean  Hs  from  HARBD,  barge  basin  east  corner,  Station  8 
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Figure  G4.  Mean  A7sfrom  HARBD,  barge  basin  north,  Station  27 


Figure  G5.  Mean  Hsfrom  HARBD,  barge  basin  south  corner,  Station  7 


Appendix  G  Plan  4c  Impact  on  Barge  Basin  Wave  Climate 


G5 


Figure  G6.  Mean  /-/sfrom  HARBD,  barge  basin  east  corner,  Station  8 


Figure  G7,  Mean  Hsfrom  HARBD,  barge  basin  north,  Station  27 
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both  plans,  with  Plan  4c  showing  a  small  improvement  over  Plan  la.  The  more 
exposed  barge  north  location  failed  to  meet  the  criterion  in  both  plans. 


Table  Gl 

Mean  and  Maximum  H&  Based  on  Numerical  Modeling 

Location 

Station 

Mean  H# ft 

Maximum  Ha,  ft 

Plan  la 

Plan  4c 

Plan  la 

Plan  4c 

Channel,  500  ft  seaward  of  jetty  tip 

17 

1.57 

1.77 

5.8 

6.6 

Channel,  beginning  of  confined  entrance 

9 

1.28 

1.31 

6.2 

5.4 

Barge  basin,  north 

27 

0.69 

0.69 

3.4 

3.5 

Barge  basin,  east  corner 

8 

0.46 

0.46 

2.3 

2.1 

Barge  basin,  south  corner 

7 

0.39 

0.39 

2.3 

1.9 

Table  G2 

Percent  Occurrence  of  Hs  Greater  than  0.3  m  (1  ft),  Based  on 
Numerical  Modeling 

Location 

Percent  Ha  >  0.3  m  (1ft) 

Plan  la 

Plan  4c 

Barge  basin  north,  Station  27 

16.5 

16.5 

Barge  basin  east  corner,  Station  8 

6.6 

5.8 

Barge  basin  south  corner,  Station  7 

5.0 

3.5 

Storm  Waves  in  Barge  Basin 

Physical  model  experiments 

A  total  of  44  wave  cases  from  the  previous  physical  model  study1  were 
analyzed  to  assess  the  impact  of  Plan  4c  relative  to  Plan  la  on  storm  waves  in  the 
barge  basin.  Physical  model  wave  gauges  in  the  barge  basin  were  positioned  in 
the  north  part  and  in  the  south  comer  (Stations  27  and  7  in  Figure  Gl).  In 
Plan  la,  Station  27  is  directly  exposed  to  waves  from  241-  to  276-deg  azimuth, 
an  arc  of  35  deg  that  includes  waves  from  the  west.  In  Plan  4c,  waves 
approaching  more  northwesterly  than  260-deg  azimuth  are  partially  blocked  by 
the  jetty. 


1  Briggs,  M.,  Lillycrop,  L.,  Harkins,  G.,  Thompson,  E.,  and  Green,  D.  (1994).  “Physical  and 
numerical  model  studies  of  Barbers  Point  Harbor,”  Technical  Report  CERC-94-14,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 
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The  44  wave  conditions  analyzed  included  eight  measured  field  conditions 
and  36  empirical  wave  cases.  Selection  criteria  for  the  eight  field  cases  ( BPS 
series  in  Briggs  et  al.  1994 *)  included  the  following: 

a.  Obtaining  the  largest  wave  heights  and  a  representative  range  of  wave 
period  and  direction  within  model  constraints. 

b.  Preference  given  to  time  after  marina  opened  in  July  1989  and  second  SX7 
directional  gauge  installed. 

c.  Maximum  number  of  operational  field  gauges  for  comparisons. 

The  eight  field  cases  had  peak  wave  periods  ranging  from  6  to  18  sec,  significant 
heights  ranging  from  2.1  m  (7  ft)  to  3.0  m  (0  ft),  approach  directions  of  208-  to 
266-deg  azimuth,  and  directional  spreading  up  to  19  deg.  All  eight  cases 
represent  rare  events  because  of  their  large  wave  heights.  The  water  level  for  all 
experiments  was  mean  lower  low  water  (mllw). 

The  36  empirical  wave  cases  had  periods  of  7.7, 11.1,  and  16.7  sec, 
significant  height  of  7.6  m  (25  ft),  and  wave  directions  of  210-,  235-,  and 
270-deg  azimuth.  These  cases  were  representative  of  very  rare  storm  events 
because  of  their  large  wave  heights.  Half  of  the  wave  cases  were  directionally 
spread  with  a  relatively  wide  spread  of  10  deg  (BPD  series),  and  half  were 
unidirectional  with  a  narrow  spread  of  1  deg  (BPU  series).  Both  narrow-  and 
wide-frequency  spreading  (i.e.,  y=  3.3  and  7.0),  representative  of  sea  and  swell, 
respectively,  were  also  included.  The  water  level  included  a  storm  surge  and 
tide  of  4.5  ft. 

Thus,  the  44  cases  for  physical  modeling  represent  moderate  to  extreme  storm 
wave  heights  and  a  representative  range  of  wave  periods,  directions,  and 
directional  spreading.  The  eight  field  cases  had  significantly  wider  directional 
spread  than  any  of  the  36  empirical  cases.  Physical  model  wave  heights  are 
based  on  data  sampling  of  each  wave  case  for  a  prototype  time  of  5.7  hr  (40  min 
model  time).  Reported  wave  heights  are  significant  heights  calculated  from 
measured  surface  elevation  time  series. 


Relative  impacts  of  Plan  4c 

Relative  performance  between  Plan  4c  and  Plan  la  was  assessed  by 
comparing  measured  Hs  values  at  the  two  barge  basin  wave  gauges  for  each  of 
the  44  wave  cases.  A  value  less  than  1.0  (or  100  percent)  indicates  that  Plan  4c 
gives  reduced  wind-wave  energy  in  the  barge  basin,  a  desirable  effect.  Overall 
barge  basin  response  is  presented  as  an  average  of  the  two  barge  basin  gauges 


1  Briggs,  M.,  Lillycrop,  L.,  Harkins,  G.,  Thompson,  E.,  and  Green,  D.  (1994).  “Physical  and 
numerical  model  studies  of  Barbers  Point  Harbor,”  Technical  Report  CERC-94-14,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 
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(numerical  model  station  numbers  are  used  to  conveniently  designate  physical 
model  gauge  locations): 


Barge  basin  response 


0.5 


* 


Hs  Station  27,  Plan  4c 

H  Station  27,  Plan  la 
v  5 


II  Station  7,  Plan  4c 

+  — - 

Hs  Station  7,  Plan  la^ 


The  overall  basin  response  indicates  a  reduction  of  storm  wave  height  in 
Plan  4c  relative  to  Plan  la  (Table  G3).  The  impact  of  Plan  4c  is  small  for  the 
empirical  cases  but  large  for  the  field  cases,  possibly  because  of  the  large 
directional  spread  used  in  the  field  case  simulations. 


Table  G3 

Average  Barge  Basin  Hs  Response  to  Storms,  Based  on  Physical 

Model  Experiments 

Experiment  Series 

Response1 

(percent) 

Field  cases 

76 

Directional  empirical  cases 

93 

Unidirectional  empirical  cases 

95 

1  Average  for  Plan  4c  divided  by  average  for  Plan  la,  converted  to  percent. 

For  the  field  experimental  series,  the  lowest  and  highest/^  response  ratios 
were  65  and  82  percent  (Figure  G8,  and  Table  G4).  The  lowest  case  approached 
parallel  to  the  entrance  channel  with  relatively  narrow  directional  spreading. 

The  highest  case  approached  the  barge  basin  more  directly  and  had  wider 
directional  spreading.  The  second  lowest  response  was  for  a  case  with  wave 
period  of  only  7.7  sec  and  approach  direction  of  218-deg  azimuth,  oblique  to  the 
barge  basin. 

For  the  directional  and  unidirectional  empirical  cases,  the  lowest  and  highest 
Hs  response  ratios  were  recorded  by  the  same  two  cases.  The  lowest  values  were 
73  and  82  percent  for  the  directional  and  unidirectional  cases,  respectively 
(Figures  G9  and  G10  and  Table  G4).  The  highest  values  were  109  and 
107  percent  for  the  directional  and  unidirectional  cases,  respectively.  The 
lowest  case  was  for  a  relatively  short  wave  period,  narrow  frequency  spreading 
(i.e.,  y  =  7.0),  and  southerly  wave  direction  (not  directly  into  the  barge  basin). 
The  highest  case  was  for  relatively  long  wave  period,  wide  frequency  spread 
(i.e.,  y  =  3.3),  and  westerly  wave  direction  (directly  into  the  barge  basin).  It 
appears  that  the  barge  basin  response  increases  as  wave  period,  frequency 
spreading,  and  barge  basin  exposure  increase. 

Differences  between  the  field  and  empirical  case  groups  can  be  attributed  to 
several  factors.  Field  cases  are  a  small  number  of  events  that  represent  a  limited 
sampling  of  conditions.  Directional  spreading  in  most  field  cases  was 
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Figure  G8.  Physical  model  Hs  ratios  for  field  cases;  directional  spreads  are  narrow 
(7-9  deg)  and  wide  14-19) 


Table  G4 

Physical  Model  Cases  with  Lowest  and  Highest  Barge  Basin 
Response 

Case 

(Plan  4c  HJ  / 

(Plan  la  percent 

Ha 

TP 

sec 

8,  deg 
az. 

Dir. 

spread 

deg 

Freq. 
spread,  y 

Field  Cases 

Lowest 

65 

2.2  m  (7.1  ft) 

14.2 

225 

9 

Highest 

82 

2.3  m  (7.4  ft) 

10.0 

243 

16 

Empirical  Cases 

Lowest: 

directional 

unidirectional 

73 

82 

7.6  m  (25  ft) 

7.6  m  (25  ft) 

7.7 

7.7 

210 

210 

10 

1 

77 

Highest: 

directional 

unidirectional 

109 

107 

7.6  m  (25  ft) 

7.6  m  (25  ft) 

16.7 

16.7 

270 

270 

10 

1 

3.3 

3.3 

significantly  greater  than  even  the  directional  empirical  cases.  A  disproportion¬ 
ate  number  of  field  cases  have  wave  approach  directions  nearly  parallel  to  the 
channel  (five  out  of  eight  field  cases  are  within  10  deg  of  being  parallel  with  the 
channel). 

Physical  model  results  also  provide  insight  on  variation  of  storm  wave 
conditions  across  the  barge  basin  (between  the  northeast  and  southwest  areas). 
Differences  between  the  two  areas  were  investigated  by  computing  the  ratio  of 
Hs  for  Plan  la  and  Plan  4c  at  Station  27  and  this  ratio  of  Hs  at  Station  7.  The 
average  ratio  for  field  cases  is  1.71  in  Plan  la  and  1.77  in  Plan  4c,  indicating  that 
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re  G9.  Physical  model  Hs  ratios  for  empirical  cases;  effect  of  frequency  and 
directional  spreading,  averaged  over  all  approach  directions 
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Figure  G1 0.  Physical  model  Hs  ratios  for  empirical  cases;  effect  of  wave  approach 
direction  and  directional  spreading,  averaged  over  all  frequency 
spreadings 
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Station  27  tends  to  have  significantly  higher  wave  conditions  than  Station  7  for 
these  storm  events.  In  empirical  studies  representing  very  rare  storms,  average 
ratios  for  directional  cases  are  1.48  for  Plan  la  and  1.47  for  Plan  4c  and  for 
unidirectional  cases  are  1.33  for  Plan  la  and  1.27  for  Plan  4c.  Thus,  wave 
heights  for  extreme  storm  events  become  less  variable  over  the  barge  basin. 
Considering  that  numerical  model  results  showed  a  ratio  of  1.77  between  mean 
Hs  values  at  Station  27  and  Station  7  (Table  Gl),  there  is  a  clear  tendency  for 
wave-height  differences  across  the  barge  basin  to  decrease  when  more  extreme 
events  “flood”  the  area  with  wave  energy. 

The  deeper  channel  in  Plan  4c  tends  to  allow  more  wave  energy  to  pass 
through  the  entrance  to  Barbers  Point  Harbor.  However,  the  jetty  helps  to  shelter 
the  entrance,  especially  for  waves  approaching  from  the  west  or  north  of  west. 
The  net  effect  on  mean  significant  wave  height  in  the  barge  basin  is  no  change 
from  the  existing  harbor.  Plan  4c  provides  a  small  improvement  in  protecting 
the  barge  basin  from  most  storm  waves. 
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Appendix  H 

Plans  7a  and  7b  Harbor 
Oscillations1 


Harbor  oscillation  studies  of  Barbers  Point  Harbor  with  a  planned  183-m- 
(600-  ft-)  wide  by  335-m-  (1,100-ft-)  long  expansion  were  requested  by  the 
Pacific  Ocean  Division  to  document  harbor  response  and  ensure  that  the  plan 
harbor  is  not  prone  to  excessive  long-period  resonance.  Two  depth 
configurations  were  studied,  based  on  physical  model  experiment  results 
(Table  HI). 


Table  HI 

Barbers  Point  Harbor  Plans  7a  and  7b 

Plan 

Entrance  Channel 

Deep-Draft  Harbor  Depth 

Width 

Depth 

7a 

140  m  (450  ft) 

13.4  m  (44  ft) 

12.8  m  (42  ft) 

7b 

140  m  (450  ft) 

14.3  m  (47  ft) 

13.7  m  (45  ft) 

Approach 

Harbor  oscillation  studies  were  performed  with  the  numerical  model 
HARBD.  The  model  was  calibrated  to  field  measurements  in  previous  studies  of 
Barbers  Point  Harbor  (Briggs  et  al.  1994).  Also,  the  model  was  used  previously 
to  study  oscillations  in  other  plan  configurations  for  the  harbor.  The  intent  in 
this  study  was  to  use  the  calibrated  model  to  study  the  new  Plans  7a  and  7b 
configurations.  The  U.S.  Army  Engineer  Waterways  Experiment  Station  has 
gained  considerable  experience  and  understanding  of  HARBD  since  the  previous 
studies  were  performed.  Experience  has  been  obtained  from  both  basic  experi¬ 
ments  with  idealized  harbors  and  practical  experience  with  other  harbor  studies 
accompanied  by  field  measurements,  including  Kaumalapau  Harbor,  Lanai,  and 


1  Appendix  H  was  prepared  by  Dr.  Edward  F.  Thompson  and  Mr.  Doyle  L.  Jones. 
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Kahului  Harbor,  Maui  (Smith  1998;  Thompson  et  al.  1996).  Based  on  these 
experiences,  it  was  evident  that  some  parameter  settings  used  in  the  earlier 
Barbers  Point  Harbor  numerical  model  study  should  be  modified  for  best  results. 

The  numerical  model  grid  for  the  existing  Barbers  Point  Harbor  developed  by 
Briggs  et  al.  (1994)  was  used  in  this  study  to  recalibrate  with  revised  parameter 
settings.  Field  data  described  by  Briggs  et  al.  (1994)  and  Lillycrop  et  al.  (1993) 
were  again  used  for  calibration.  Parameters  for  the  recalibrated  numerical  model 
are  given  in  Table  H2.  Key  differences  in  this  study  include  a  constant  bottom 
friction  value  and  a  harbor  boundary  reflection  coefficient  less  than  1.0. 


Table  H2 

Parameter  Values  Used  in  HARBD 

Parameter 

Value 

Bottom  friction,  /? 

0.032 

Coastline  reflection,  Krcoast 

1.0 

Depth  in  infinite  region,  dfar 

8.5  m  (28.0  ft) 

Harbor  boundary  reflection,  Kr(wave  periods  of  25-100  sec) 

0.95 

Harbor  boundary  reflection,  Kr( wave  periods  >  100  sec) 

0.96 

Finite  element  numerical  grids  representing  Plans  7a  and  7b  were  developed 
by  modifying  the  Plan  4c  grid  used  previously.  Modifications  include  reducing 
the  size  of  the  harbor  expansion  to  183  m  (600  ft)  by  335  m  (1,100  ft),  moving 
the  channel/harbor  depth  transition  to  the  area  where  the  channel  opens  into  the 
harbor  basin  (as  recommended  from  physical  model  studies),  and  modifying 
project  depths  as  in  Table  HI.  Plan  7a  and  7b  grids  are  similar  except  for  project 
depths  imposed  in  the  harbor  and  entrance  channel.  Grid  characteristics  are 
summarized  in  Table  H3. 


Table  H3 

Grid  Size 

Harbor  Plan 

Number  of: 

Elements 

Nodes 

Solid  Boundary  Nodes 

Semicircle  Boundary  Nodes 

7a,  7b 

12,501 

6,603 

592 

112 

Wave  periods  ranging  from  25  to  1,000  sec  were  modeled  using  parameter 
settings  derived  by  recalibration  (Table  H2).  Incident  long  wave  height  was 
0.06  m  (0.2  ft),  representing  a  realistic  long  wave  event,  and  approach  direction 
was  approximately  perpendicular  to  shore.  Harbor  response  information  was 
saved  at  32  station  locations  (Figure  HI). 
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Figure  HI .  Output  station  locations 


Results 

Harbor  oscillation  results  are  summarized  as  amplification  factor  versus 
wave-frequency  plots  for  each  station.  Amplification  factor  is  defined  as  long 
wave  amplitude  at  the  station  location  divided  by  incident  long  wave  amplitude. 
Plots  for  representative  stations  in  the  deep-draft  harbor,  barge  basin,  and 
channel  are  included  in  this  report  (Figures  H2-H15).  Plans  7a  and  7b,  included 
together  in  the  plots,  had  very  similar  responses.  Resonant  peaks  in  Plan  7b 
were  shifted  to  slightly  higher  frequencies  relative  to  Plan  7a  because  of  the 
deeper  project  depths. 

Except  for  very  low  frequencies  at  which  the  entire  harbor  rises  and  falls  in 
unison  (frequencies  less  than  0.002  Hz;  periods  longer  than  500  sec), 
amplification  factors  in  the  channel  and  main  harbor  (Stations  1-6  and  28)  are 
nearly  all  less  than  2.0,  with  a  maximum  of  3.0.  In  the  deep-draft  harbor 
expansion  area  (Stations  29-32),  a  tendency  for  stronger  amplifications  is 
evident.  Comer  Stations  30  and  31  are  most  active,  though  maximum 
amplification  factors  are  still  around  3.0  or  less.  Barge  basin  comers  (Stations  7 
and  8)  are  more  active  than  the  deep-draft  harbor.  Amplification  factors  exceed 
3.0  for  several  peaks,  but  are  always  less  than  4.0. 
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Figure  H4.  Wave  amplification  factors,  Station  3 


Figure  H5.  Wave  amplification  factors,  Station  4 
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Figure  H8.  Wave  amplification  factors,  Station  7 
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Figure  H9.  Wave  amplification  factors,  Station  8 
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Figure  HI 2.  Wave  amplification  factors,  Station  29 
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Figure  HI 3.  Wave  amplification  factors,  Station  30 
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Past  experience  with  physical  and  numerical  modeling  of  harbor  oscillations 
suggests  the  following: 

a.  If  the  amplification  factor  near  a  pier  at  any  relevant  long  wave  period  is 
greater  than  5,  some  operational  difficulties  may  be  encountered. 

b.  If  the  amplification  factor  is  greater  than  10,  major  operational 
difficulties  can  be  expected. 

By  these  criteria,  the  channel,  the  deep-draft  harbor,  and  the  barge  basin  in 
Plans  7a  and  7b  are  not  prone  to  excessive  long-period  resonance.  A  recent 
analysis  of  long  wave  velocities  in  Barbers  Point  Harbor  with  335-m  (1,100-ft) 
by  335-m  (1,100-ft)  expansion  (Plan  4c)  shows  a  high-energy  long  wave  event 
may  produce  velocities  of  over  0.5  m/sec  (1  knot)  adjacent  to  the  expansion 
comer  where  it  opens  into  the  main  harbor  basin  (Thompson,  Boc,  and  Nunes 
1998).  These  localized  velocities  may  be  a  concern  in  Plans  7a  and  7b. 
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d.  The  transition  should  be  moved  to  the  opening  of  the  harbor  basin.  At  the  proposed  transition  location,  there  is 
less  wave  energy;  if  vessel  shear  occurs,  the  harbor  pilot  would  have  more  room  to  react/correct.  Model  navigation 
study  data  support  a  0.6-m  (2-ft)  transition. 

Based  on  numerical  investigations,  the  deeper  channel  in  Plan  4c  tends  to  allow  more  wave  energy  to  pass 
through  the  entrance;  however,  the  jetty  helps  to  shelter  the  entrance,  especially  for  waves  approaching  from  the  west 
or  north  of  west.  The  net  effect  on  mean  significant  wave  height  in  the  barge  basin  is  no  change  from  the  existing 
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deep-draft  harbor,  and  barge  basin  in  Plans  7a  and  7b  are  not  prone  to  excessive  long-period  resonance.  Another 
recent  analysis  of  long-wave  velocities  for  Plan  4c  shows  a  high-energy  long-wave  event  may  produce  velocities  of 
over  0.5  m  (1  knot)  adjacent  to  the  expansion  corner  where  it  opens  into  the  main  harbor.  These  localized  velocities 
may  be  a  concern  in  Plans  7a  and  7b. 


